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ABSTRACT 


This  study  Is  intended  to  develop  a  new  concept  in  modeling  a  mon- 
orgonized  system  and  to  initiate  implementation  and  testing  of  this  model.  It 
was  desired  to  not  only  have  the  capability  to  predict  the  output  of  the  organiza¬ 
tion  but  to  also  be  able  to  characterize  the  internal  operation  of  the  system. 

These  objectives  may  be  accomplished  by  combining  the  best  features  of  modeling 
from  an  economist's  viewpoint  and  a  control  system  engineering  standpoint.  The 
economists  consider  primarily  the  input  and  output  flows  and  the  control  system 
theory  dwells  on  the  consideration  of  parameters  influencing  system  performance. 

] 

The  use  of  the  state  variable  equations  permitted  the  accomplishment  of  these  i 

I 

objectives. 

The  state  variable  equations  were  rewritten  in  a  difference  form  and  com¬ 
bined  into  a  matrix  form.  Discrete  data  from  the  real-world  system  is  used  to 

generate  a  matrix  which  is  characteristic  of  the  system.  This  fitting  procedure  ^ 

is  a  unique  feature  of  the  Dynamic  Organizational  Network  Analysis  (DONA) 

method.  Multiple  regression  analysis  is  used  in  the  DONA  method  to  develop 

the  system  description  matrix.  This  matrix  is  the  primary  component  of  the 

DONA  model.  The  DONA  model  will  produce  system  outputs  and  future  system 

state  variables  given  the  system  inputs. 

Since  the  DONA  method  of  producing  a  model  is  o  new  approach,  a 
question  arises  as  to  its  validity.  The  usual  procedure  to  test  ihe  validity  of  a 
model  is  to  use  historical  data  from  the  real-world  and  compare  the  model  outputs 

ii 


with  the  results  of  the  simulated  organization.  This  method  does  not  permit  the 
desired  control  over  the  stimulus  data.  A  laboratory  concept  was  conceived 
to  validate  the  DONA  mudel  and  to  discover  its  basic  characteristics;  i.e. 
frequency  response,  stabdity,  etc.  The  standard  for  comparison  was  a  computer 
simulation  of  a  sales  company.  This  simulation  contained  known  operational 
capabilities  and  was  therefore  a  good  model  to  use  for  the  comparison. 

The  complete  DONA  method,  DONA  mode!  ond  laboratory  concept  were 
written  Irrto  u  computer  program  using  FORTRAN  language. 

The  results  of  the  program  pointed  out  some  unsuspected  problems  in 
modeling  a  mun-organized  system.  The  computation  of  the  system  characteristic 
matrix  leads  to  a  better  understanding  of  the  system  operation.  The  computer 
program  also  contains  a  unique  feature  of  printing  out  how  well  the  real-world 
system  has  been  described  in  the  DONA  model. 

Since  this  study  was  the  initial  effort  in  developing  the  DONA  concept, 
there  were  several  areas  discovered  thot  needed  further  study.  The  thesis  con- 
eludes  with  some  recommendations  for  further  exploration. 
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CHAPTER  I 


GENERAL  MODEUNG 

Introduction 

During  the  post  decade,  significant  advancements  in  business  organiza¬ 
tion  modeling  have  been  achieved.  This  may  be  attritutabic  to  increased 

1* 

interest  in  systems  analysis  and  to  the  advent  of  the  high-speed  computer. 

The  concept  of  modeling  is  based  on  an  obstraction,  sirr^lification, 

or  idealization  of  the  system  or  event.  The  model  may  be  used  for  a  variety 

of  purposes.  The  most  common  of  these  is  prediction  or  optimization  of  the 

system  being  modeled.  Depending  on  th«  approach  used  in  constructing  the 

model,  another  extremely  valuable  use  of  the  model  Is  to  describe  the  octivity 
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going  on  within  the  system.  R.  A.  Brown  in  "Modeling  the  Orgonizotion" 
discusses  modeling  the  orgonizotion  as  seen  from  two  viewpoints-economists 
and  management.  The  first  viewpoint  is  that  of  an  accountant  where  the 
system  is  represented  as  consisting  of  inputs  and  outputs.  From  the  economist's 
opproach,  very  little  is  said  about  what  occurs  internol  to  the  system.  The 
second  viewpoint  is  that  of  the  control  system  engineer.  The  system  is  viewed 
from  within.  Control  of  the  organization  is  attempted  through  the  manipulotion 

*  Numbers  in  superscript  throughout  the  thesis  indicate  the  reference  as  listed 
in  the  List  of  References. 
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of  selected  internal  characteristic  quantities  to  have  the  greatest  effect  on 
the  system  performance.  Then  the  control  systems  approach  to  modeling  accom¬ 
plishes  the  desired  goals  of  modeling,  predicting  future  events  and  optimizing 
system  performance,  and  further  permits  an  understanding  of  the  endogenous 
(internal)  variables  that  control  the  operation  of  the  organization. 

The  Problem 

With  the  understanding  of  why  the  modeling  of  an  organization  is 
performed,  it  is  pertinent  to  direct  attention  to  one  type  of  an  organization — 
the  man-organized  system.  The  man-organized  system  is  characterized  by 
elements  such  as  people,  material,  money,  and  information.  The  set  of  elements 
together  with  the  relationships,  define  the  man-orgonized  system. 

The  man-organized  system  to  be  considered  in  this  thesis  is  a  lorge 
governmental  agency  In  general,  the  method  of  modeling  developed  herein 
is  not  restricted  to  this  type  of  organization.  It  will  be  seen  that  the  method 
may  be  applicable  to  other  organizations. 

In  the  modeling  of  this  large  governmental  agency,  it  was  required 
that  the  modeling  effort  be  performed  in  a  reasonobly  short  time.  The  model 
must  characterize  the  system  through  use  of  elemental  quantities  related  to 
basic  system  operation.  Every  variable  in  the  system  need  not  be  accounted 
for  in  the  model.  This  results  in  a  gross  representation  of  the  system  and  this 
was  the  desired  approach. 

Before  pursuing  a  new  approach  to  the  problem,  it  is  informative  to 
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note  how  the  problem  is  currer\tly  being  handled.  The  objectives  ore  to  know 
something  about,  (1)  how  the  system  will  behave  in  the  future  given  that 
estimates  are  possible  of  the  workloads  applied  to  the  organization,  (2)  what 
parameters  in  the  system  have  the  greatest  control  over  the  system  performance 
and  (3)  what  is  the  magnitude  of  system  variation  with  a  known  change  in  the 
significant  porameters.  The  usual  procedure  is  to  collect  information  from 
selected  monitoring  points  in  the  system  and  assemble  this  informotion  into 
a  report,  usually  a  recurring  report.  Using  the  data  from  these  reports,  an 
analysis  is  made  and  an  experience  factor  is  applied  to  estimate  the  internal 
system  characteristics  and  the  system  outputs. 

Another  method  used  in  attacking  the  problem  and  searching  for  an 
answer  to  the  objectives  is  to  model  the  organization.  Modeling  permits  the 
representation  of  the  system  in  quantifiable  terms.  However,  the  complexity 
of  a  large  organization  mokes  it  very  difficult  to  represent  the  system  in  mathe¬ 
matical  form.  Further,  the  solution  of  the  mathemoticol  expressions  may  be 
even  more  elusive.  The  concept  of  computer  simulation  has  been  used  in  solving 

this  problem.  One  of  the  most  extensive  works  on  computer  simulation  of  large 
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organizat ions  is  by  Jay  W.  Forrester.  Forrester  introduces  a  different  concept 
in  modeling  large  organizations.  He  uses  information  feedback  ideas  to  describe 
"industrial  dynamics,"  his  concept  of  dyfiomic  modeling.  Forrester  olso  develops 
his  own  computer  longuage,  known  as  DYNAMO,  to  use  in  his  simulations.  The 
Industrial  Dynamics  method  of  modeling  is  suited  to  on  effort  which  requires  de¬ 
tailed  representation  of  a  system.  This  method  may  require  years  for  accomplish- 


ment  of  fhe  objectives. 

Since  it  wcs  desirable  to  improve  upon  the  current  methods  of  collecting 
and  analyzing  reports  or  using  the  lengthly  Industriol  Dynamics  simulation 
procedure,  o  search  began  for  a  new  method  of  modeling  that  would  be  timely 
and  lend  to  the  understanding  of  the  performance  of  the  man-organized  system. 

This  thesis  will  give  insight  to  a  new  method  of  modeling  and  develop 
the  concept  into  the  early  stages  of  computer  program  implementation  and 
testing.  Chapter  II  gives  on  example  of  a  computer  simulation  of  a  large  system. 
This  simulation  shows  that  mathematical  expressions  are  required  to  use  the 
simulation  method.  This  simulation  was  included  to  illustrate  the  complexity 
of  this  approach  even  for  a  small  organization.  Chapter  111  introduces  the  new 
concept  for  modeling.  The  Dynamic  Orgonizotion  Network  Analysis  (DONA) 
method  is  developed  in  this  chapter.  Chopter  IV  discusses  o  unique  procedure 
in  validating  a  model.  The  laboratory  concept  for  the  testing  of  a  system  model 
is  introduced.  Chapter  V  gives  the  results  of  this  study  and  provides  areas 
for  further  study  in  the  use  of  the  DONA  method. 


CHAPTER  I] 


MAN-ORGANIZED  SYSTEM  AS  SIMULATED  WITH  A 
COMPUTER  MODEL 

Design  and  Development  Considerations  in  Simulation 

This  chapter  describes  a  computer  model  which  is  used  to  simulate  a 
man-organized  system.  This  model  serves  two  purposes.  It  serves  to  illustrate 
the  detail  and  logic  needed  to  construct  a  computer  model  and  it  will  be  used 
in  the  validation  procedure  described  in  further  chopters. 

"Modeling"  and  "simulation"  are  often  times  used  interchangeably.  There 
is  a  difference  in  these  terms.  The  model  as  used  herein  means  the  set  of 
mathematical  equations  which  quantitatively  describes  the  system  being  modeled. 
A  computer  model  is  the  mathematicol  model  written  in  a  form  to  be  used  by 
the  computer.  Then  simulation  is  the  use  of  the  model.  The  methodology  of 
simulation  begins  w ith  constructing  the  mathematical  model.  In  describing  a 
large  governmental  agency  of  a  thousand  or  more  people  in  a  quantitative  form 
the  task  may  become  quite  formidable.  In  fact  the  expression  of  certain  actions 
in  mathematical  equations  moy  be  elusive  or  impracticol.  The  detail  required 
for  such  a  characterization  will  be  illustrated  in  the  model  described  below. 

The  model  will  assume  certain  stochastic  processes.  Although  the  model  is 
hypothetical,  these  stotistical  quantities  could  be  verified  if  the  model  were 
really  describing  a  real  system;  i.e.,  the  probability  distribution  and  its  mean 
and  variance  of  a  particular  vorioble  could  be  derived  from  the  real  system. 
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Another  Qssunrption  in  the  hypothetical  model  is  that  the  policies  of  the  or¬ 
ganization  ore  correctly  defined  mothemotically.  Of  course,  one  use  of  the 
model  is  to  test  the  optimality  of  the  policies.  The  assumptions  and  policies 
ore  used  to  bound  the  problem.  This  model  is  therefore  an  effort  to  reduce  a 
complex  system  to  manageable  proportions. 

Description  of  SIMCO 

The  computer  model  is  a  modified  model  that  was  originally  developed 
4 

by  McMillan  and  Gonzalez.  This  model  is  known  as  SIMCO  and  is  written 
in  FORTRAN  language.  It  is  a  model  for  a  total  system  with  specified  parometers. 
The  total  system  is  a  business  firm  which  sells  one  product.  It  does  not  manu¬ 
facture  the  product  but  does  order  it  from  the  supplier.  In  this  sense  SIMCO 
is  a  sales  distributor  and  is  characterized  in  general  as  o  warehousing  situation. 
Customers  place  demands  in  a  random  manner  and  their  orders  may  be  filled 
immediately  or  not  filled  depending  on  avoilobil ity  of  stock  .  The  suppi  ier 
requires  lead  time  for  delivery  which  is  also  o  random  process.  The  business 
activity  is  simulated  on  a  weekly  basis.  The  model  performs  colculotlons 
of  the  activities  at  the  end  of  each  week  ond  stock  is  received  at  the  end  of 
eoch  week. 

The  computer  progrom  describes  the  random  processes  and  company 
policies  of  SIMCO  first.  These  processes  and  policies  are  described  as  follows. 


Lead  time  for  delivery  from  the  supplier  varies  from  four  weeks  to  eight 
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weeks  with  the  following  probability  distribution: 

Lead  Time  Probability 

(weeks) _ 


4  .10 

5  .15 

6  .50 

7  .15 

8  .10 

Customer  demand  is  uncertain  from  week  to  week  but  is  normally  distributed 
with  a  mean  of  50  units  per  week  and  a  standard  deviation  of  15  units. 

SlMCO  management  allows  a  price  markdown  if  weekly  sales  have 
not  been  deemed  satisfactory.  The  normal  retail  price  is  $100  per  unit; 
however,  if  demand  during  the  week  was  less  than  70  per  cent  of  forecasted 
demand,  the  retail  price  is  marked  down  to  $85.  The  price  morkdown  is 
also  dependent  on  another  condition.  The  inventory  on  hand  at  the  end  of 
the  week  plus  new  stock  scheduled  for  delivery  over  the  weekend  must  be 
more  than  150  per  cent  of  the  forecasted  demand  for  the  coming  week.  This 
"overstock"  condition  will  cause  price  markdown.  Forecasted  demand  is 
found  by  averaging  the  damand  during  the  past  three  weeks.  Experience  shows 
that  the  price  markdown  has  caused  a  sales  stimulation  of  5  per  cent. 

At  the  end  of  the  week  SlMCO  determines  the  activity  of  the  past 
weeks  and  then  places  its  order  for  new  stock.  The  amount  to  order  is  deter¬ 
mined  by  several  fc.o'  >rs,  one  of  which  is  the  desired  inventory  level.  The 
desired  inventory  is  three  times  the  demand  forecast  for  the  coming  week.  The 
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amount  to  order  is  the  sum  of  the  quantity  of  units  sold  during  the  post  week, 
one  half  of  the  difference  between  the  desirea  and  actual  inventory,  and  one 
half  of  the  difference  between  the  desired  and  actual  supply  backorder.  The 
desired  supply  backorder  is  the  e>(pected  demand  during  the  lead  time,  based 
on  the  demand  forecast  for  the  coming  week. 

SlMCO's  experience  has  shown  that  about  20  per  cent  of  tl^e  backorders 
ore  lost  if  they  are  not  filled  the  following  week.  Thus  a  backorder  situation 
is  permitted  in  SIMCO. 

The  wholesale  cost  for  units  remains  constant  at  S50  per  unit.  The 
payment  is  made  immediately  after  receipt  of  the  new  stock. 

McMillan  and  Gonzalez  continue  their  description  of  SIMCO  by 
next  looking  at  the  money  transactions.  There  are  two  types  of  costs,  fixed 
and  variab'le.  Fixed  costs  are  $2,000  per  week.  Variable  costs  are  20  per 
cent  of  the  investments  in  inventory  at  the  beginning  of  the  week,  at  a  yearly 
rate,  and  8  per  cent  of  the  investment  in  accounts  receivable  at  the  beginning 
of  the  week,  at  a  yearly  rate.  The  fixed  and  variable  costs  are  cosh  costs 
and  represent  an  outflow  of  cash  at  the  end  of  each  week.  Cash  inputs  to 
SIMCO  occur  with  the  sale  of  the  units.  Again  it  has  been  observed  that  one 
half  of  the  sales  ore  in  cash  and  the  other  half  are  on  credit.  The  charge  soles 
are  payable  four  weeks  later  and  do  not  incur  a  service  charge. 
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SIMCO  Sequence  of  Operations 

The  sequence  of  operaHons  for  simulating  100  weeks  operations  is  shown 
in  the  FIc  w  Diagram,  Figure  2.1.  There  is  a  moin  program  with  four  subroutines. 
Initiolization  occurs  in  the  main  program  and  the  subroutine  supply  is  c^  Med 
to  determine  how  much  new  stock  wos  received  over  the  weekend.  The  demand 
subroutine  is  then  colled.  This  segment  simulates  the  demand  that  SIMCO  will 
experience  in  the  forthcoming  week.  As  McMillan  and  Gortzalez  hove  constructed 
SIMCO,  all  cosh  transactions  occur  next  and  these  occur  In  the  moin  program. 

The  computer  model  concludes  with  o  subroutine  to  forecast  demand  and  a  sub¬ 
routine  to  determine  the  quantity  of  units  to  be  ordered.  The  entire  program 
is  pot  in  a  loop  which  repeats  100  times,  each  loop  representing  one  week  of 
operation. 


SIMCO  Initial  Coiiditions  and  Programs  Elements 

The  programs  elements  os  summarized  above  will  now  be  described  in 
more  detail.  Before  doing  this,  the  initial  conditions  of  SIMCO  ore  considered. 
Any  time  that  o  simulation  is  run  on  the  computer,  there  is  the  problem  of  starting 
conditions  and  equilibrium.^  If  the  starting  conditions  vary  greatly  from  the 
steady  state  conditions,  the  time  to  reach  stabilization  may  be  quite  long. 

]t  is  important  to  reach  steady  state  conditions  because  the  information  content 
of  the  simulation  i:,  contained  In  the  equilibrium  state.  The  transient  state 
nioy  be  greatly  shortened  by  the  choice  of  reasonable  starting  conditions.  The 
following  initial  conditions  were  used  in  the  SIMCO  model. 
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Demand  during  the  past  10  weeks  has  been 
50  units  each  week; 

SIMCO  is  destined  to  receive,  on  eoch  of  the 
coming  six  weekends,  50  units  of  new  stock; 

The  actual  inventory  on  hand  is  150  units; 

Demand  forecast  for  the  current  week  is  50  units; 

There  are  no  customer  backorders  outstanding; 

There  is  no  price  markdown  planned  for  the 
current  week; 

$2,500  is  scheduled  to  be  received  from  accounts 
receivable  collections  over  each  of  the  coming  four 
weeks  (thus  total  accounts  receivable  ore  $10,000); 

Cash  on  hand  stands  at  $4,000. 

It  was  stated  that  the  Supply  subroutine  was  called  to  determine  the 
amount  of  stock  received  over  the  weekend.  Since  orders  have  already 
been  placed,  SIMCO  new  stock  may  be  "destined  to  be  delivered"  in  1 ,  2,  3, 
4,  5,  6,  7,  or  8  weeks.  These  simulated  deliveries  in  the  pipeline  will  be 
designated  by  the  variable,  PIPLIN(K).  A  rondom  number  generator  will 
simulate  the  lead  time  according  to  the  probability  distribution  listed  above. 
Subroutine  Random  produces  two  random  numbers.  The  random  numbers,  RND2, 
is  uniformly  distributed  from  0  to  1  and  is  used  in  the  Supply  subroutine.  The 
random  number  DEVlOT,  is  normally  standard  distributed  (mean  =0,  Standard 
deviation  =  1  .)  and  is  used  in  the  subroutine  Demand.  The  quantity  of  units 
ordered,  DO,  of  the  Order  subroutine  is  applied  to  the  lead  time.  In  effect. 
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the  Supply  subroutine  creates  o  delay,  an  uncertain  delay,  In  receiving  the  units 
ordered.  The  output  of  the  Supply  subroutine  is  the  stock  received,  SR.  The 
main  program  takes  the  stock  received  and  fills  the  customer  backorders,  if  any. 
An  inventory  updating  occurs  by  adding  to  the  actual  inventory  at  the  end 
of  the  previous  week,  AIEW,  the  stock  received  over  the  weekend  to  yield 
the  actual  inventory  beginning  the  new  week,  AIBW.  The  customer  backorder 
queue,  CBQ,  is  subtracted  to  yie  d  the  adjusted  actual  inventory  beginning 
the  new  week,  AAIBW.  At  this  point,  units  of  lost  sales,  ULS,  may  occur  if 
the  customer  backorders  exceed  the  inventory,  a  stockout  condition.  This  con¬ 
dition  is  noted  in  the  program. 

The  main  program  calls  the  subroutine  Demand  to  determine  the  demand 
which  SlMCO  will  experience  during  the  coming  week.  A  rondom  number 
generator  produces  the  random  numbers  according  to  a  normal  distribution  and 
the  demand  is  simulated.  Using  the  price  markdown,  if  any,  which  was  deter¬ 
mined  in  the  main  program  through  consideration  of  the  management  policy 
for  the  price  markdown,  the  subroutine  Demand  applies  the  sales  stimulation 
factor,  SSF,  to  increase  the  simulated  demand  by  5  per  cent. 

The  main  program  then  begins  a  series  of  update  transactions.  Using 
the  current  week's  demand,  ADDW,  the  inventory  on  hand  after  customer 
backorders  were  filled  at  the  beginning  of  the  week,  AAIBW,  and  backorder 
sales  at  the  beginning  of  the  week,  BS,  the  following  are  calculated: 

Total  Sales  during  the  week,  TUS 

Actual  Inventory  at  the  end  of  the  week,  AlEW 
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Gross  Profit  and  Gross  Income  from  the  week's 
soles,  WGP  and  GIS,  respectively. 

The  main  program  continues  with  the  money  transactions.  The  Cosh  on  Hand, 
COH,  will  be  increased  by  accounts  receivable  collected  during  the  week  and 
income  for  cosh  sales  during  the  week.  COH  is  decreased  by  payment  for 
new  stock,  fixed  charges  of  $2000  per  week,  and  variable  costs.  The  variable 
chorges  consist  of  the  Inventory  bolding  charge  and  the  investment  charge  on 
the  total  accounts  receivable.  The  accounts  receivable  for  a  particular  week 
is  one-half  the  gross  income  from  sales  four  weeks  past.  The  total  accounts 
receivable  is  the  sum  of  accounts  receivable  from  sales  1,  2,  3,  and  4  weeks 
past.  The  weekly  net  profit,  WNP,  is  then  determined  from  the  weekly  gross 
profit  minus  the  operating  costs  (fixed  costs  and  variable  costs). 

The  demand  forecasting  subroutine  merely  takes  the  demand  during  the 
most  recent  past  three  weeks  and  averoges  them  to  produce  the  forecasted 
demand  for  the  ccming  week.  This  process  of  determining  the  forecasted  demand 
can  be  recognized  os  o  moving-average,  smoothing  operation. 

The  main  program  uses  the  actual  demand  during  the  week  and  the 
forecasted  demand  to  implement  the  management's  policy  of  declaring  or 
not  declaring  a  price  markdown. 

The  Reorder  subroutine  determines  the  number  of  units  to  be  ordered 
from  the  supplier  or  the  end  of  the  week.  Knowing  that  the  desired  inventory 
is  three  times  the  forecasted  demand,  the  difference  of  actual  inventory  ond 
desired  inventory  is  determined  os  the  first  element  of  the  decision  to  reorder. 
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The  second  element  is  based  on  the  difference  of  the  desired  and  actual  supply 
backorder  pipeline.  The  quantity  of  units  to  order  is  then  the  total  units 
sold  during  the  week  plus  one-half  the  sum  of  the  two  differences  noted  above. 

The  SlMCO  total  systems  simulation  is  concluded  with  a  report  gen¬ 
erating  sequence.  McMillan  and  Gonzalez  recognized  the  need  to  produce 
something  that  could  be  used  to  analyze  the  changes  In  total  operation  when 
testing  the  value  of  a  single  management  policy.  To  this  end  the  main  program 
contains  a  series  of  statements  to  determine  averages  and  variances  for  all 
significant  data.  These  are  printed  out  for  analysis. 

The  complete  SIMCO  computer  listing  is  shown  in  Appendix  A.  The 
output  from  this  simulation  is  also  shown  in  Appendix  A. 

Modificotion  to  SIMCO 

It  is  readily  apparent  that  the  SIMCO  model  described  above  is  extremely 
rich  in  its  formulation  and  operation.  The  model  contoins  o  v/eolth  of  features 
related  to  mathematical  operotions.  These  are  queues,  logs,  smoothing  averages, 
and  forecasting  built  into  the  model.  However,  the  use  that  this  model  is  to 
be  put  to  later  in  this  project  requires  some  additional  features.  SIMCO  is 
a  single  product  model.  Since  it  is  a  single  produc*  model,  it  obviously  has 
no  interactions  with  other  products  and  this  is  desired.  SIMCO  does  not 
contain  personnel  considerations.  It  is  desired  to  simulate  personnel  through 
consideration  of  skills  and  numbers  of  personnel .  With  this  in  mind,  extensive 
modifications  were  made  to  SIMCO  as  shown  in  the  Flow  Diagram  in  Figure  2.2. 
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Figure  7.7  (Coofinued).  Tiow  Diagiam  of  Modified  SlMCO 
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A  second  product  line  was  added  to  the  SlMCO  model.  All  four 
subroutines  of  Supply,  Demand,  Forecasting,  and  Reorder  hod  to  be  developed 
for  the  second  product.  Now  customers  con  place  demand  for  product  one 
ond  product  two,  or  both.  The  demand  chorocteristics  may  be  the  so.  ie  or 
difVrrent.  The  modification  will  permit  independent  operations  of  the  two  pro¬ 
ducts. 

A  completely  new  subroutine  was  developed  for  personnel  considerations. 
A  hypothetical  staffing  and  manning  organization  was  used  to  bose  the  skills 
and  numbers  of  personnel.  Part  of  the  original  SlMCO  fixed  costs  were  allo¬ 
cated  to  salaries  for  the  SlMCO  manager  ond  administrative  personnel.  These 
costs  and  personnel  remained  fixed  regardless  of  company  operotions.  However, 
additional  personnel  were  added  as  salesmen.  The  salesmen  were  assumed  to 
have  a  certain  productivity;  i.e.  the  salesmen  hod  the  copobility  of  selling 
up  to  20  units  per  week  of  their  particular  product.  The  number  of  solesmen 
for  each  product  was  based  of  the  forecasted  demand  for  the  coming  week. 

This  resulted  in  a  hire-fire  operat  ion.  But  to  moke  the  model  o  little  more 
realistic,  the  fire  procedure  caused  SlMCO  to  pay  two  weeks  solory  for 
severance  pay. 

The  personnel  subroutine  provided  an  excellent  opportunity  to  bring  in 
some  interaction  between  the  product  lines.  The  personnel  subroutine  wos 
further  developed  to  allow  transfer  of  solesmen  from  one  product  line  to  the 
other.  This  permits  o  salesman  to  be  transferred  from  a  product  whose  demand 
is  forecasted  low  to  a  product  whose  demond  is  forecasted  high.  Personnel 
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ocHons  now  include  hire-fire-transfer  actions.  This  policy  has  desirable  side 
effects.  Rather  than  having  such  high  variations  in  hire-fire,  the  added  transfer 
policy  has  dampened  these  variations  and  has  resulted  in  o  more  realistic 
company  simulation. 

Modified  SIMCO  Results 

The  modified  SIMCO  may  be  used  in  the  same  manner  as  the  original 
SIMCO.  The  modified  SIMCO  may  be  used  to  test  various  management  policies. 
Since  it  was  not  desired  to  use  the  modified  SIMCO  for  this  purpose,  this 
was  not  pursued.  The  flow  diagram  and  the  complete  computer  listing  for  the 
modified  SIMCO  is  shown  in  Appendix  B.  It  is  sufficient  to  note  that  the  modi¬ 
fied  SIMCO  did  produce  results  to  verify  that  the  modifications  were  effective 
in  that  a  second  product  line  could  be  simulated  independent  of  the  first  product 
line  except  for  the  desired  interaction  of  personnel  considerations.  The  personnel 
modification  further  expanded  the  usefulness  of  the  simulation. 


HAPTER  III 


DYNAMIC  ORGANIZATIONAL  NETWORK  ANALYSIS 
Introduction 

This  chapter  will  describe  the  Dynamic  Orgonizational  Network  Analysis 
(DONA)  method  and  will  describe  the  use  of  this  method  as  applied  to  modeling 
the  man-organized  system, 

''he  DONA  method  is  a  new  approach  to  the  modeling  of  socio-economic 

systems.  The  model  is  derived  from  control  theory  considerations  much  the  same 

os  used  by  the  electrical  engineers  for  analysis  of  systems.  The  strength  of  this 

approach  lies  in  the  fact  that  the  system  is  considered  from  both  the  economic 

and  management  view  point.  The  method  for  producing  the  DONA  model  adopts 

the  regression  technique  from  the  economic  models  ond  preserves  the  multi- 

2 

dimensionality  of  the  management  models. 

The  control  theory  approach  will  be  exploited  in  the  manner  described 

6  7  8  9 

by  H.  E.  Koenig  '  '  '  when  he  used  the  state  variable  equations  to  model  a 

'0 

socio-economic  system,  specifically  a  university.  Koenig,  et.al.,  in  his 

8 

text  ''Analysis  ot  Discrete  Physical  Systems  "  stated,  "Machine  simulations 
os  well  as  recent  advances  in  the  theory  of  stability,  control,  optimization, 
design,  and  synthesis  are  all  based  on  what  are  referred  to  os  stote-space  models." 
The  state-space  model  concept  is  the  basis  of  the  DONA  model.  The  approach 
applies  to  linear  systems  but  this  limitation  does  not  invalidate  its  application  to 
non-linear  systems. 

This  thesis  will  extend  the  extensive  work  by  Koenig  by  applying  the 
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regression  technique  to  discover  the  "structure"  of  the  system  to  be  modeled. 

Koenig  derived  the  structure  of  the  university'*^  through  the  lengthly  process 

ofwriting  the  mathematical  equations  for  each  act ivity  w ithin  'he  university. 

He  odmits  that  all  variables  were  not  quantifiable.  It  will  be  shown  that  the 

DONA  model  does  not  require  the  representation  of  the  system  as  a  set  of 

mathematical  equatio'is.  The  structure  of  the  system  can  be  determined  from 

observations  by  regression  methods.  The  odvantcges  of  the  DONA  model  will 

lie  in  the  ease  of  estimation  of  the  required  constants,  representing  the  structure, 

because  of  the  existence  of  an  objective  method  for  determining  those  constants. 

The  resulting  representation  will  be  meaningful  because  it  can  be  reloted  to 

2 

control  theory,  which  will  give  guidance  for  system  improvements. 

The  Stote-Space  Approach 

8  1 1 

In  using  the  state-space  approach  '  to  systems  ancly..is,  it  is  necessary 
to  be  able  to  represent  the  system  with  equations.  Such  a  representation  is  the 
mothemotical  model  and  forms  the  starting  point  for  the  problem  of  system  analysis. 
Aiilock  diagram  representation  for  o  system  is  shown  in  Figure3. 1  . 


X(t)  . 

System 

Y{t)  . 

V(t),  Hi) 

Figure  3 . 1 


Block  Diagram  of  a  System 
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In  Figure  3  .1  the  X(t)  is  o  vector  representing  the  inputs  or  stimuli! 
and  the  Y(t)  is  a  vectoi  representing  the  outputs  or  responses.  The  c(t)  and 
i  (t)  are  vectors  representing  the  state  variables  and  the  first  derivative  of  the 
state  variables,  respectively.  Then  the  specification  of  the  values  of  the  state 
variables  at  some  time  instoni  is  the  state  of  the  system  at  that  Instant. 

Assuming  a  general  n-th  order  linear  system,  the  mother  atical  model 
is  formulated  from  a  set  of  simultaneous  first-order  differential  equations  and  a 
set  of  simultaneous  algebraic  equations.  Differential  equations  to  represent  the 
general  'r+h  order,  linear,  time-varying  system  for  a  single  input  and  output, 
can  be  v^ritten  as 

dV*) 

dt 

d  vjit) 

-  =p2iW  V|(t)  +P22(0  ■‘■P2n^*^  'n^^^ 

dt 

d  i  (t) 

-  ^  +  •••  +  P„.(0  ^(0  +q„(0  ^(0- 

nil  n^  z  nn  n  n 

dt 

The  output  equation  is  in  algebraic  form 

y(t)  =  mj(t)  ;j(t)  +  m2(t)  ';f2(0  +  •  •  •  '■^(t)  +  n(t)  x(t)  (3-2^ 

In  equations  (3-1)  and  (3-2)  the  pjj(t),  q.(t),  m.(t),  ond  n(t)  are  t  ime-vorying 
coefficients  of  their  respective  variables. 

The  above  equations  considered  only  o  single  input  and  single  output. 
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However,  it  is  possible  to  represent  a  system  with  several  inputs  and  outputs 
by  o  set  of  equations  by  adding  to  each  equation  the  appropriate  contribution 
from  each  input.  Similarly,  the  output  equation  is  expressed  In  terms  of  the 
state  variables  by  a  sot  of  equations  instead  of  a  single  equation. 

For  a  large  complex  system,there  is  a  need  for  a  compact  notation  when 
there  are  many  inputs,  outputs,  and  state  variables.  Through  the  use  of  vectors 
and  matrices,  the  equations  may  be  reduced  to  a  compact  form. 

If  there  arc  n  state  variables,  these  con  be  combined  in  a  column 
vector,  or  an  (n  x  1)  matrix,  and  designated  os  ^  (t).  Thus, 

(n  X  1)  (3-3) 

m  input  signals,  the  input  vector 

(m  X  1)  (3-4) 

Likewise,  if  there  are  p  output  signals,  the  output  vector  is 


Ttf)  = 


i,(0 


>(0 


is  known  as  the  state  vector.  If  there  are 
will  be  designated  as 


X(t)  - 


’‘|(0 


^2(0 


X  (0 
rn 
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y,(0 

Y{0  =  I  I  (p  1)  (3-5) 


yp(t) 


This  thesis  will  consider  only  the  time-invariant  system  which  results  in  constant 
coefficients  in  the  equations.  The  coefficients  of  the  state  variables  can  be 


represented  by  an  (nxn)  matrix. 
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The  matrix  to  represent  the  input  coefficients 

is 

■^11 

q,2 

^Im 

Q  = 

‘’21 

<^22 

'^2m 

'’nl 

•  •  •  •  •  • 

^nm 

(n  X  n) 


(n  X  m) 


(3-6) 


(3-7) 


For  the  output  equations,  the  matrix  of  constant  coefficients  for  the 


state  variables  is 
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M  = 
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"'l2 
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(p  X  n) 


(3-8) 


The  final  matrix  represents theinput  coefficients  which  are  a  part  of 
the  output  equations. 
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(3-9) 


By  substituting  the  above  defined  vectors  and  matrices  into  equations 
(3“1)  and  (3-2)  the  following  equations  result  and  are  known  as  the  state  variable 
equations . 


i  {»)  =p  -no  +Q  x(t) 

Y(t)  v(t)  +N  X(t) 


or 


'rjt)' 

’  P 

Q  ' 

’  Mt)' 

Y(t)_ 

M 

N 

_X(t)_ 

(3-10) 

(3-11) 

(3-12) 


The  use  of  the  state  variable  equations  will  require  the  equations  to 
utilize  discrete  information.  This  will  require  a  difference  form  of  the  state 
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voriablo  equations.  The  continuous  form  of  the  equations  (3-12)  are  converted 
to  the  difference  form  as  follows 


=  p  4(t)  +QX(t) 
h 

dL  V(l)]_  lim  4(t  +h)  -  liiftl 

^  Q - - W 

dt  h 

V  (t  +  h)  =  hP  'i(t)  +  V  (t)  +  hQ  X(t) 


Then 


Define  the  S  matrix, 

f]'?-’  si 


(3-13) 


(3-U) 


(3-15) 


(3-16) 


(3-17) 


and  substituting 


i(t  +  1) 
_  y(0_ 


(3-18) 


This  equation  is  in  a  form  that  suggests  the  use  of  multiple  regression 
analysis  (see  equation  (10)  in  Appendix  C).  In  fact  the  S  matrix  may  be  deter¬ 
mined  if  regression  methods  ore  applied  to  the  system  to  be  modeled.  If  equation 
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(3-18)  were  used  to  model  the  system,  it  is  suspected  that  this  would  predict 

system  performance  in  both  short  term  ond  long  term  time  durations.  These 

choracteristics  may  be  desirable  in  some  modeling  problems  but  this  study  is 

more  concerned  with  short  period  predictions,  especially  because  of  the  assumed 

constant-coefficient  form  o.  the  equations.  Another  form  of  the  state  equations 

appears  to  be  applicable  to  the  short  term  prediction.  Further  study  is  needed 

to  understand  the  characteristics  of  using  the  model  in  equation  (3-18)  and  the 

short  term  model  to  be  described  below.  It  will  be  based  on  differences  of  a 

variable  from  one  time  period  to  another.  This  is  equivalent  to  high-pass  filter- 

12 

Ing  of  the  data . 

The  term  "forward  difference"  refers  to  the  variable's  value  in  one  time 
period  in  advonce  minus  the  variable's  curie.nt  value.  In  mathemotical  form 
the  forward  difference,  ^(x),  is 

‘l‘(x)  =x(t  +  1)  -  x(t).  (3-19) 

Similarly,  the  backword  difference  is  the  variable's  current  value  minus  its  value 
one  time  perirxd  in  the  past.  The  backward  difference, ^(x) ,  is 

v(x)  ^  x(t)  -  x(t  -  1)  .  (3-20) 

Monipulotion  and  substitution  of  equation  (3-  and  (3-20)  into  equation 

(3-16)  yields 


(3-21) 
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The  methodology  for  developing  the  DONA  model  is  succinctly 
summarized  inequations  (3-19),  (3-20),  and  (3-21).  The  forward  ond  backward 
differences  are  determined  for  the  appropriate  variables  selected  from  the  system 
to  be  modeled.  Linear  algebra  is  used  to  solve  the  matrix  equation  (3-21)  for  the 
S  matrix  and  this  is  the  matrix  which  characterizes  the  system. 

The  DONA  model  is  represented  by  equation  (3-18).  Having  determined 
the  S  matrix  as  described  above,  the  DONA  model  is  determined. 


CHAPTER  IV 


THE  LABORATORY  EXPERIMENTAL  SYSTEM 

Validafing  The  Model 

In  the  previous  chapters,  two  completely  unrelated  types  of  models  for 
man-organized  systems  were  shown.  A  new  modeling  procedure  based  on  the 
state-space  approach  was  discussed  ond  a  computer  simulation  of  an  organization 
was  explained.  This  chapter  will  show  o  way  to  validate  a  model  which  is  quite 
different  from  the  "classical"  method  of  validating  a  model.  It  will  be  seen 
that  this  new  procedure  provides  irrany  important  advantages  over  the  usual  model 
volidatlon  procedure. 

Validation  or  verificotion,  as  it  sometimes  colled,  ^  is  a  procedure 
used  to  determine  the  "truth"  or  "worth"  of  a  model.  Then  the  concept  for 
volidotion  is  based  on  establishing  o  set  of  criteria  for  differentiating  between 
what  is  true  and  whot  is  not  true.  This  is  not  unlike  theproblem  of  testing  a 
statistical  hypolhesis.  However,  the  procedure  or  the  method  foi  validation  is 
the  question.  The  most  common  method,  the  "classical"  method,  of  validating 
the  model  is  to  use  data  from  the  real  world,which  the  model  is  simulating,  and 
observe  the  performance  of  the  model.  As  pointed  out  in  Chapter  1,  a  model 
may  serve  two  purposes.  It  moy  be  used  to  determine  optimal  policies  and 
operation  of  a  system  or  it  may  be  used  to  predict  the  performance  of  the  system. 
Naylor,  et.  al.  ,  '  in  discussing  the  problem  of  verification  said: 
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It  is  our  position  that  tne  ultimate  test  of  a  computer 
simulation  model  is  the  degree  of  accuracy  with  which  the  model 
predicts  the  behavior  of  the  actual  system  (which  is  being  sim¬ 
ulated)  in  the  future. 

This  stotement  is  true  but  it  does  not  fulfill  the  requirement  for  validating 
the  DONA  model.  An  immediate  requirement  exists  for  testing  the  DONA 
model  and  the  method  for  producing  the  DONA  model  to  discover  some  of  the 
basic  characteristics.  This  requires  control  over  the  variables  being  applied  to 
the  model.  The  actual  system  cannot  be  manipulated  to  produce  the  range  and 
magnitude  of  variables  needed  to  discover  the  operating  characteristics  of  the 
DONA  model.  For  this  reason,  a  laboratory  concept  was  developed  to  validate 
the  method.  In  this  "laboratory,"  the  parameteis  can  be  controlled  to  determine 

the  range  and  responsiveness  of  the  model.  The  laboratory  concept  was  used  and 
will  be  described  in  the  following  sections. 

The  Laboratory  System 

The  concept  for  validation  of  the  DONA  Model  and  hence  the  quality 
of  the  valid  method  used  to  produce  the  DONA  model  was  to  use  SIMCO  to 
simulate  the  "real  world"  organization  and  implement  the  DONA  methodology 
to  model  SIMCO.  SIMCO  was  used  because  its  characteristics  were  known. 
Using  the  same  inputs  to  SIMCO  and  DONA  the  validity  of  DONA  con  be  esti¬ 
mated  by  comparing  the  outputs  of  the  two  models.  The  complete  "laboratory" 
system  is  shown  in  Figure  4. 1 .  The  implementation  of  this  laborotory  was  through 
the  use  of  a  computer  program  in  FOf-'RAN  V  language  to  describe  the  elements 


aborotory  System 
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and  flov/s  in  the  laboratory.  The  computer-based  laboratory  system  in  Figure 
4.2  will  now  be  described  in  detail. 

SIMCO  was  described  in  detail  in  Chapter  II  but  some  additional  con¬ 
siderations  had  to  be  employed  to  use  the  simulation  in  the  laboratory  system. 
blMCO  produces  about  50  variables.  Each  of  these  variables  had  to  be 
analyzed  to  group  them  in  one  of  three  types  of  variables — input,  output,  or 
state  variables.  Even  though  the  SIMCO  simulation  ,  per  se,  is  unconcerned 
about  which  variable  is  which,  the  identification  of  these  variables  for  the 
state-space  approach  is  very  important. 

It  is  not  sufficient  to  draw  an  imaginary  line  around  SIMCO  and 
identify  everything  crossing  the  line  as  either  an  output  or  input.  This  is  the 
accountant's  viewpoint  of  identifying  the  flows  as  inputs  or  outputs.  The 
control  theory  viewpoint  must  also  be  e;^loited.  The  system  must  still  be 
defined.  The  same  imaginary  line  may  be  used,  but  now  the  terminals  must 
be  considered  as  either  a  propensity  variable  or  flow  rate  variable  depending 
on  its  characteristics.  This  analysis  performed  on  SIMCO  yielded  some  rather 
unexpected  results.  Even  though  there  were  many  variables  in  SIMCO,  the 
terminal  analysis  showed  that  there  were  not  enough  variables  produced  by 
the  simulation.  Some  equations  were  then  introduced  in  the  simulation  to 
complete  the  terminal  requirements.  An  eximple  of  this  problem  is  the 
quantity  of  stock  received  each  week,  a  flow  rote  variable,  from  the  supplier. 
The  stock  received  is  an  input  to  SIMCO;  however,  there  is  no  singularly 
identifiable  propensity  at  the  output  which  is  complementary  to  this  input. 

A  new  output  variable,  wholesale  cost,  was  developed  to  account  for  the  com- 
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plementary  input  of  stock  received.  In  this  cose,  the  wholesale  cost  was  $50 
per  unit  and  so  the  total  wholesale  cost  wos  equol  to  the  stock  received  timts 
$50.  This  equation  was  then  put  in  the  program.  Moving  completed  this  input- 
output  analysis,  the  remaining  v  ariables  were  exomined  for  their  appropriateness 
of  describing  the  state  of  SIMCO.  The  remaining  voriobles  were  subsequently 
identified  and  labeled  as  '*ate  variables.  For  the  purposes  of  testing  the  labora¬ 
tory  system,  all  variables  were  not  used.  The  variables  selected  for  the  checkout 
were  2  input,  2  output,  and  5  state  variables.  The  input  variables  were  the 
actuol  demands  during  the  week  for  product  I  and  for  product  2,  (ADDW  and 
ADDWl).  The  state  variables  were  units  ordered  for  products  1  and  2  (LiO  and 
U01),  actual  inventory  end  of  week  for  product  1  (AIEW),  and  forecasted  demand 
for  product  1  and  2  (FD  and  FDl).  The  output  variables  were  salaries  for  per¬ 
sonnel  on  product  1  (SAL)  and  wholesale  cost  of  product  1  received  from  supplier 
(WC). 

The  analysis  of  the  "real  world"  model  was  not  completed  with  the 
identifying  of  the  input,  output  and  state  variables.  The  DONA  methodology 
requires  a  matrix  inversion  at  a  point  in  its  sequence  of  operations.  The  elements 
in  the  matrix  to  be  inverted  were  initially  derived  from  o  series  of  values  from 
the  "real  world"  inputs,  outputs  and  state  variables.  A  motrix  inversion  requires 
that  the  matrix  be  non-singu'ar.  Specificolly,  the  matrix  must  not  contain  any 
rows  that  are  linear  combinations  of  onother  ond  there  must  not  be  a  column  or 
row  with  all  zeroes  for  the  values.  The  computer  program  for  the  loborato'y 
showed  that  SlMCO  was  indeed  producing  both  of  these  objectionable  traits 
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and  thereby  causing  niatrix  singularity.  This  says  that  SIMCO  in  simulating 
the  real  world  was  producing  redundant  information.  Linear  combinations  of 
some  variables  were  appearing  in  the  SIMCO  simulation.  The  other  untenable 
situation,  all  zeroes  in  a  row  or  column,  resulted  from  the  characteristics  of  the 
DONA  methodology.  The  method  to  produce  the  DONA  model  requires  that 
the  variables  not  remain  at  a  constant  value  or  that  the  slope  (first  derivative) 
not  be  constant. 

The  point  to  be  made  in  using  SIMCO  or  using  a  real  world  organization 
is  that  a  careful  analysis  must  be  performed  in  selecting  the  variables  and  noting 
the  values  taken  by  the  variables.  The  computer  program  has  several  built-in 
features  to  detect  the  matrix  singularii  ies  but  the  laoeling  of  variables  as  input, 
output  or  state  must  be  a  manual  task. 

To  continue  the  description  of  the  laboratory  system,  the  DONA 
methodology  block  in  Figure  4. 1  contains  the  multiple  regression  analysis  procedure 
for  developing  the  "S"  matrix.  In  the  laboratory  system  the  "S"  matrix  will  char¬ 
acterize  the  SIMCO  Soles  Company.  The  computer  program  for  this  segment  of 
the  laboratory  system  has  some  unique  features  some  of  which  were  dictated  by 
the  limitations  of  the  computer  (UNIVAC  1 108)  and  some  added  to  produce 
traceability  ond  analysis  of  the  computations  used  in  the  DONA  methodology. 

These  features  may  be  recognized  by  referring  to  the  computer  flow  diagram 
in  Figure  4.3.  The  computer  program  for  this  segment  of  the  laboratory  system 
is  called  DONAFT  (DONA  Fit)  and  implements  the  DONA  methodology. 

The  DONAFT  subroutine  begins  by  odmitting  the  preselected  input, 
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output  and  state  variables  of  SIMCO  generated  for  82  weeks  of  SIMCO  opera¬ 
tions.  The  appropriate  forward  and  bockword  differences  were  calculated  for 
all  observations  (82  weeks).  The  forward  ond  backward  differences  as  defined 
in  Chapter  111  were  the  difference  in  value  between  the  present  time  and  one 
time  past  for  the  backward  difference  and  the  difference  between  one  time 
interval  in  the  future  and  the  current  time  for  the  forwotd  difference.  The 
backward  differences  of  the  input,  output  and  state  variables  were  determined 
and  the  forward  differences  of  only  the  state  variables  were  determined.  These 
differences  were  performed  ot  each  time  interval  beginning  with  week  2  and 
ending  in  week  81 .  Weeks  1  and  8'*  had  to  be  available  because  in  the  series 
of  differences,  the  backward  difference  requires  one  interval  in  the  post,  week 

1  when  week  2  Is  considered.  Similarly,  the  forward  difference  requires  one 
intervol  in  advance  of  the  series,  week  82  when  week  81  is  considered.  The 
results  of  the  differences  govt  j  sets  of  differences  for  each  of  the  2  input, 

2  output,  and  5  state  variables. 

The  array  of  differences  must  be  formed  for  the  next  process  in  o  manner 
such  that  the  top  half  rows  of  the  motrix  contain  the  independent  variables. 

Each  column  of  the  matrix  represents  a  time  interval.  This  matrix  is  called 
the  data  mati  ix  or  "X"  matrix. 

The  "X"  matrix  is  transferred  into  the  ViULTR  (Multiple  Regression) 
subr'outine.  The  mean  of  each  variable  was  determined  and  subtracted  out 
leaving  the  deviation  from  the  mean  for  each  variable.  The  deviations  from 


the  mean  are  still  in  the  difference  form.  The  deviations  do  not  reflect  back 
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to  the  real  value  of  the  individual  varioble.  All  computations  are  based  on 
the  forward  and  backward  differences.  Having  the  deviations  from  the  mean, 
the  sum  squares  are  computed  for  all  combinations  of  the  variables.  The  "A" 
matrix  and  "K"  matrix  as  described  in  Appendix  C,  have  now  been  determined. 

The  computer  program  produces  these  sum  squares  in  a  square  array 
AUK  (A  union  K).  In  actuality,  the  square  A  matrix  is  the  upper  left  quarter 
of  the  square  AUK  matrix  and  the  square  K  rtrotrix  is  the  upper  right  quarter  of 
the  AUK  matrix.  The  usual  procedure  would  be  to  take  the  A  matrix  into  a 
matrix  inversion  subroutine  to  perform  the  inversion.  However,  one  of  the  afore¬ 
mentioned  problems  now  arise.  Assuming  the  A  matrix  has  hod  the  zeroes 
removed,  if  they  we^e  present  in  a  whole  row  or  column,  and  assuming  linear 
combinations  are  not  present  in  the  A  matrix,  the  nature  of  the  A  motrix  couses 
the  inversion  process  to  exceed  the  computer  capability.  Specificolly,  the  A 
matrix  contains  values  ranging  from  10^  to  10  ^  when  using  the  SlMCO  date. 
Even  if  the  program  were  written  in  double  precision,  it  is  doubtful  that  a  good 
inversion  could  be  obtained.  The  multiplicotion  processes  in  the  motrix  inversion 
would  exceed  the  capability  of  the  computer  registers.  The  answer  to  this 
problem  was  to  "normalize"  the  AUK  matrix.  This  con  legitimately  be  done 
through  conversion  of  the  A  matrix  into  a  matrix  of  intercorrelations.  The 

elements  of  the  AUK  matrix  represent  nO..:  i.e.  the  covarionce  times  the  number 

'J 

of  sample  points. 

An  element  of  the  AUK  matrix  equals  the  sum  of  products  (x.  •  x.). 

Then  the  covariance  results  from 


39 


n 

V 

k  =  1 


(x.|  X.. )  =  n^.. 

ik  jk  ij 


(origin  at  the  means) 


(4-1) 


wh_.c:  J..  is  defined  as 
■J 


'j 


k  -  1 


n 


n 
I, 
k  = 


n 


(4-2) 


The  ith  element  on  the  main  diagonal  (ith  column  and  ith  row)  of  the 


AUK  matrix  is  nO,.,  also  written  nC."  ,  The  same  is  true  of  the  jth  element;  i.e. 
nOjj  =  .  The  AUK  matrix  has  been  written  into  the  AAA  matrix  so  that  the 


AAA  motrix  may  be  used  as  a  working  matrix. 

The  intercorrelctions  or  Peorson  product -moment 
are  defined  as  ^ 
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This  cquotiori  rriOy  be  used  by  modification  os  fcilows 
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By  applying  this  relatior.  tothe  AAA  matrix  the  resulting  matrix,  the 
intercoi  re  lot  ion  matrix,  ranges  in  value  from  -I  to  +1.  The  intercorrelation 
matrix  was  then  applied  to  o  matrix  inversion  subroutine  wttich  ciso  had  on 
Inlerncl  matrix  inversion  occurocy  improvement  routine.  Tlie  inverse  is  returned 
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to  the  form  of  sum  of  products  through  an  inverse  transformation  using  the  inter- 
correlation  relation  (4-4)  in  reverse.  The  A  inverse  matrix,  submatrix  of  the 
inverted  AUK  matrix,  is  pre-r.iultiplied  with  the  K  matrix  in  conformance  with 
the  multivariate  regression  analysis  (sec  Appendix  C).  The  result  is  the  CUB 
(G  onion  B)  matrix  where  the  B  matrix,  sub-matrix  of  CUB,  is  the  matrix  of 
coefficients  of  regression  (Appendix  C). 

This  concludes  the  operation  of  MULTR  subroutine  except  for  a  very 
helpful  set  of  statements  which  compute  and  printout  the  statistical  properties 
for  each  dependent  variable  as  a  function  of  the  independent  variables.  One 
important  printout  is  the  amount  of  variability  that  has  been  explained  in  the 
dependent  varioble  through  the  independent  voriobles. 

The  DONA  methodology  and  the  DONAFT  subroutine  is  completed  with 
one  more  process.  The  "S"  matrix  as  defined  in  Chapter  111,  equation  (3-17  ) 

is  obtained  by  taking  the  right  half  of  the  CUB  matrix,  the  C  sub-motrix,  end 
transposing  it.  The  resulting  S  mat'ix  characterizes  SlMCO  Sales  Company, 

The  DONA  model  block  in  Figure  4-1  represents  the  DO^'IA  model  as 
discussed  in  Chapter  111.  The  jro;ess  in  this  s»-age  is  depicted  by  the  following 
matrix  equations 


;  (r  +  I  ) 
_Y(t)  _ 


(4-5) 


The  left  side  of  the  equotion  is  the  output  from  the  DONA  model. 

The  TONA  subroutine  of  the  laboratory  system  computer  program  (see 
Figure  4.4)  performs  the  calculations  of  the  DONA  model  using  equation  (4-5). 
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Figure  <.4  FlO'V  Diogram  of  DON  A  Subroutine 
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The  S  matrix  was  developed  in  the  DONAFT  subroutine  and  entered  in  the 
DONA  subroutine.  The  state  variable  initial  conditions  used  for  SIMCO  are 
also  imposed  on  the  DONA  subroutine  for  the  first  time  interval  only.  After 
the  first  time  interval  the  'i'(t  +  1),  produced  as  an  output  of  the  DONA  model, 
are  used  as  the  input  state  variables  for  the  next  time  interval.  The  X(t) 
inputs  are  provided  the  DONA  model  for  each  time  interval.  The  outputs 
of  the  DONA  model  are  the  state  variable  vector  for  the  future  time  interval 
and  the  output  vector  for  the  simulated  system. 

The  outputs  of  SIMCO  and  DONA  are  finally  compared  as  shown  in 
Figure  4-1 .  The  comparison  is  mode  on  all  parameters  desired  or  deemed  appro¬ 
priate  for  consideration. 

Operation  of  the  Laboratory  System 

The  operation  of  the  computer-based  laboratory  system  is  one  complete 
fiow;  however,  the  sequential  actions  may  be  portrayed  as  occuring  in  two  major 
parts.  Phase  1  Is  the  initial  phase  with  the  purpose  of  determining  the  S  matrix, 
the  system  characterization  matrix.  The  flow  through  the  laboratory  system  for 
this  phase  is  shown  in  Figure  4.5  by  the  solid  lines.  In  Phase  1  the  data  is 
taken  from  the  actual  system  (SIMCO  in  this  case)  and  entered  into  the  DONAFT 
subroutine  where  the  DONA  methodology  is  accomplished.  The  result  of  this 
process  is  the  S  matrix. 

Phase  II  in  the  operation  of  the  laboratory  system  is  the  performance  of 
the  DONA  mode!  and  comparison  of  the  outputs.  Having  determined  the  S  matrix 
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in  Phase  1,  the  DONA  model  is  complete  and  ready  to  function.  Phose  II  may 
be  considered  as  the  "free-running"  state  of  the  laboratory  system  in  that  both 
models  are  continuously  producing  outputs  and  state  voriables.  This  phase  is 
shown  in  Figure  4.  6  by  the  solid  lines.  It  will  be  noted  that  the  DONAFT 
subroutine  has  been  "disconnected"  in  this  phase  of  the  laboratory  operation. 

Thvj  laboratory  system  has  great  flexibility  primarily  because  it  is  a 
computer  progrcim  and  changes  moy  easily  be  made.  This  flexibility  may  be 
illustrated  in  another  way.  Suppose  that  after  having  disconnected  the  DONA 
fit  procedure  the  comparison  procedure  has  signaled  that  the  S  matrix  is  no  longer 
valid.  That  is,  the  S  matrix  no  longer  characterizes  the  system  that  it  is  repre¬ 
senting.  This  requires  generation  of  a  new  matrix  or  updating  of  the  S  matrix. 

The  laboratory  system,  perm  .  this  updating  in  a  very  uncomplicated  manner. 

The  complete  laboratory  system  is  rerun  using  new  dato  from  the  actual  system 
to  be  modeled.  The  two  phase  operation  occurs  automatically  and  thus  o  new 
system  matri,  and  its  DONA  model  are  generated. 

Validation  of  the  DONA  Method 

The  major  consideration  in  developing  and  using  the  loborotory  concept 
was  to  have  control  over  the  range  and  i/pes  of  voriobles  which  the  DONA 
model  may  experience.  With  such  exoct  control  then  some  conclusions  may  be 
reached  on  the  performance  of  the  DONA  model  and  the  method  used  to  produce 
the  model . 

The  inputs  or  stim.ulus  data  to  SIMCO  and  DONA  may  be  any  one  or  a 


Figure  4.6  Phase  II  of  the  Laboratory  System  Operation 
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combination  of  signals  composed  of  impulses,  step  functions,  sinusoids,  or 
random  noise.  In  the  early  stages  of  conducting  the  experiment,  only  one  type 
of  stimulus  was  considered.  It  was  a  normally  distributed  "noise"  with  a  mean 
of  50  ond  a  standard  deviation  of  15.  This  stimulus  was  the  customer  demand 
as  explained  in  Chapter  II  for  the  SIMCO  simulation. 

The  comparison  step  in  the  laboratory  system  was  bosed  on  computer 
printouts  of  the  SIMCO  and  DONA  state  variables  and  output  variables.  In 
a  line  by  line  comparison  of  the  SIMCO  result  and  the  DONA  result  much 
information  was  gained  on  the  performance  of  the  DONA  '  del. 

Results  of  the  Laboratory  System 

The  computer  program  for  the  laboratory  system  in  Figure  4. 1  was  written 
according  to  the  flow  diagrams  in  Figures  4.2,  4.3,  and  4.4.  The  complete 
program  listing  is  shown  in  Appendix  D.  The  program  was  checked  out  by 
running  the  individual  subroutines  ond  the  parts  tested  out  successfully  on  the 
computer  (UNIVAC  1108). 

Two  test  models  were  developed  to  discover  some  basic  behavior  char¬ 
acteristics  of  the  DONA  model.  The  first  test  model  in  Figure  4.7  contained 
only  delays.  A  discrete  approximation  of  a  sawtooth  wove  with  a  period  of  8 
time  units  was  fed  in  to  the  model.  In  one  experiment  the  sawtooth  was  delayed 
one  time  unit  followed  by  a  second  delay  of  ore  time  unit.  A  state  variable 
was  monitered  between  the  two  delays.  The  same  sawtooth  wove  was  also  fed 
in  ta  a  two  unit  time  delay.  The  above  was  repeated  but  used  an  input  of  noise 
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rather  than  the  sawtooth  wave.  The  results  were  that  a  very  high  correlation 
(100%)  could  be  obtained  by  using  one  time  unit  delay  regardless  of  input — 
noise  or  sawtooth  wave.  The  high  correlation  is  only  obtainable  by  using  a 
state  variable.  When  an  output  variable  is  used  for  correlation  after  one  time 
unit  delay,  very  poor  results  are  obtained.  An  explanation  for  the  better 
performance  of  a  state  variable  over  an  output  variable  comes  from  the  state 
equation  model.  The  state  variables  have  the  advantage  of  knowing  values 
one  time  interval  in  the  past,  whereas,  the  output  variables  are  determined 
from  current  time  calculations. 

When  a  two  unit  time  intervol  was  encountered,  poor  correlation  resulted 
with  either  input.  This  is  because  only  o  first  time  difference  is  used  in  the 
method.  Another  result  from  the  model  in  Figure  4.7  was  the  relatively  high 
correlation  (30%)  from  the  output  of  the  second  delay  when  noise  was  applied. 

In  comparison,  a  sawtooth  input  resulted  in  only  o  2%  correlation  at  the  output 
of  the  second  delay,  "^his  implies  that  the  higher  frequency  noise  is  less  filtered 
by  the  two  cascaded  deioys  than  the  lov'er  frequency  sawtooth  wave.  The  fre¬ 
quencies  of  the  noise  one  sawtooth  signals  result  from  the  first  differences  (in 
time)  taken  in  the  forward  and  bockward  differences  of  the  variables,  equations 

(3.  19)  and  (3.20). 

The  second  test  model  used  to  discover  the  basic  characteristics  of  the 
DONA  method  is  shown  in  Figure  4.8.  The  objective  was  to  examine  what 
happened  when  two  signals  were  summed.  The  result  was  tliat  there  was  100% 
correlation  of  the  output  with  the  inputs.  Further,  the  output  was  equally 
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weighted  and  weighted  to  the  maximum  extent  on  each  of  the  inputs.  This  high 
correlation  in  the  summing  operation  supports  the  theory  that  the  results  of  the 
output  variables  are  highly  significant  when  current  time  information  is  used 
rather  than  delayed  or  past  time  is  used  in  deriving  the  output  prediction. 

The  laboiotory  computer  program  was  also  successfully  operated  with 
the  SIMCO  simulation.  The  results  of  this  experiment  are  shown  in  Appendix  E. 

It  was  noted  that  the  behavior  of  the  SIMCO  predictions  was  not  good.  It 
was  expected  that  the  SIMZO  results  .vould  not  be  closely  comparable  because 
the  DONA  model  indicates  the  low  (1%  to  58%)  correlotion  of  the  variables. 

If  additional  analysis  and  adjustment  of  the  DONA  model  were  mode  to  include 
more  of  the  significant  variables,  higher  correlation  would  result.  This  was 
done  in  the  experiment  and  the  correlation  on  one  variable  was  raised  from 
23%  to  58%.  But  there  must  also  be  some  consideration  given  to  the  behavior 
of  the  DONA  as  revealed  in  the  results  of  the  test  mcdels.  Therefore,  the  con¬ 
struction  of  the  model  is  highly  significant  in  the  worth  of  its  predictive  behavior. 

Developing  the  considerations  to  be  employed  in  the  construction  of 
a  model  was  a  major  result  of  this  study.  The  DONA  model  characteristics 
dictate  that  in  modeling  the  real-world  system  the  total  delays  in  the  system 
must  be  represen: ed  in  the  DONA  model  by  cascaded  unit  deloys.  The  question 
is,  how  many  of  these  delays  ore  needed  for  the  representation?  This  is  answered 
through  the  comparison  of  the  real-world  and  DONA  values  in  the  MULTR 
Subroutine  of  the  laboratory  computer  program.  The  MULTR  subroutine  shows 
where  the  DONA  model  needs  improvement.  It  shows  which  variables  ore 
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accounting  for  the  variability  and  also  shows  ♦he  strength  of  the  variables' 
contribution . 

Another  consideration  in  model  construction  is  the  tendency  of  DONA 
to  use  "alibi"  variables;  that  is,  a  variable  completely  unrelated  to  the  operation 
or  flow  may  be  called  in  from  another  segment  of  the  model  to  explain  the  vari¬ 
ability  of  the  variable  under  study.  Again,  the  MULTR  subroutine  will  show 
the  effects  of  the  "alibi"  variables.  It  may  be  interesting  in  further  study  to 
ottempt  to  force  out  these  unwanted  variables  of  the  solution. 

Model  construction  may  be  further  enhanced  through  a  technique  to 
discover  linear  combinations.  A  problem  mentioned  earlier  in  this  chapter 
was  the  occurrence  of  linear  combinations  in  the  charocterization  of  o  man- 
organized  system.  When  a  matrix  inversion  was  attempted,  it  was  to  no  avail 
because  of  the  singular  matrix.  Linear  combinotions  in  this  connection  ore 
considered  as  the  form  where  one  equation  is  identical  to  the  other  multiplied 
by  a  constant.  It  is  recognized  that  there  are  other  forms  of  linear  combinations. 

During  the  course  of  this  research,  it  was  discovered  that  there  was  o 
way  to  recognize  the  troublesome  linear  combinations.  Consider  the  following 
matrix  as  representing  a  linear  combination  of  two  equations. 


12 

X22 

""23 

"24 

12 

0x^2 

°"24 

(4-6) 


Assume  it  is  then  placed  in  a  larger  matrix  and  made  symmetrical. 
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Linear 

combination 


"11 

''12 

0x^2 

^14 

''12 

X22 

UX22 

2 

"24 

ox, 2 

0x22 

a  X22 

°"24 

""24 

"44 

(4-7) 


Because  of  considerations  discussed  eorlier,  it  was  required  to  transform 
the  matrix  Into  intercorrelation  form.  Then  by  the  definition  of  intercorrelation, 
equation  (4-3),  the  values  in  the  main  diagonal  will  always  be  +1 .0.  The 

matrix  (4-7)  in  intercorrelation  form  will  have  a  main  diagonol  with  all  values 

2  2 
of  +1 .0;  i.e.  ,  ^  =  ^<22  ~  ^  ^^22  ~  ^^44  ~  ^  result  is  that  a  =1.0 

and  a  =+l.  In  matrix  (4-7)  a^22  equals  ^  1  since  X22  =  1-0.  The 

significance  of  this  is  that  a  linear  combination  of  the  type  described  may  be 

detected  in  the  intercorrelation  matrix  by  noting  the  presence  of  on  element 

equal  to  =  1  not  on  the  main  diagonal.  This  condition  must  exist  for  a  linear 

combination  but  it  is  not  the  only  condition.  The  inspection  for  ones  in  the 

intercorrelation  matrix  is  useful  as  a  rapid  means  to  determine  if  there  is  a 

possibility  of  linear  combinations.  The  results  of  the  laboratory  system  confirmed 

the  above.  In  Appendix  E  the  AUK  matrix,  AAA  matrix,  ond  Intercorrelation 

matrix  are  symmetrical  in  .  upper  left  quarter.  These  parts  of  the  motrices 

are  the  only  parts  used  in  the  matrix  inversion  routine.  As  expected,  the 

Intercorrelation  matrix  contains  ones  on  the  main  diogonal  and  there  are  no 

other  ones  in  the  upper  left  quarter  of  the  computer  printout.  It  will  also  be 

noticed  that  the  inverted  matrix  is  symmetrical.  In  Appendix  E  the  Improved 

Inverse  Matrix  A  (Intercorrelated  Form)  and  the  Improved  Inverse  Matrix  A 
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(No  Intercorrelations)  ore  symmetrical.  If  any  of  these  matrices  are  unsymmet- 
riccl  Joring  the  course  of  using  the  loboratory  system,  the  input  data  Is 
immediately  suspected  of  containing  irregularities. 

It  was  the  plan  of  research  to  understand  some  fundamental  characteristics 
of  the  behavior  of  the  DONA  model  and  to  understand  how  to  apply  the  method. 
This  wcs  accomplished  but  there  are  many  areas  that  need  further  exploration. 
Since  this  study  was  the  first  attempt  to  implement  the  DONA  method,  this  has 
been  only  a  beginning  in  the  research  required  to  adequately  understand  and 
apply  the  method.  It  appears  that  enough  has  been  discovered  in  this  initial 
effort  to  be  at  least  familiar  with  the  capabilities  and  limitations  in  the  use 
of  the  DONA  model. 


CHAPTER  V 


SUMMARY 

The  Approach 

This  project  resulted  from  a  need  to  model  o  large  governmentol  agency. 
The  means  of  developing  the  model  must  be  more  timely  than  previous  methods 
which  may  lequire  up  to  10  years  for  formulation.  1(  was  desired  that  the  model 
predict  in  gross  terms  as  opposed  to  detailed  prediction  of  system  variables. 

The  approach  to  this  problem  was  to  use  the  best  features  of  modeling 
from  the  economist's  view  and  the  control  theory  standpoint.  The  state  variable 
approach  was  used  to  accomplish  those  objectives.  The  use  of  state  variables 
to  produce  a  model  was  known  as  Dynamic  Organizational  Network  Analysis 
(DONA),  The  DONA  model  for  modeling  o  man-organized  system  was  developed 
from  a  methodical  means  of  evaluating  the  ri\an-organized  system.  It  is  this 
method  of  producing  the  DONA  model  that  is  the  unique  feolure  of  this  project. 

Since  the  DONA  method  of  modeling  is  a  new  approach,  a  question  arises 
concerning  Its  validity.  Therefore,  it  is  desired  to  validate  the  DONA  model 
and  this  in  turn  will  validate  the  method  used  to  produce  the  model.  A  computer- 
based  laboratory  system  was  conceived  for  validotion  purposes.  The  idea  of  a 
"laboratory"  was  used  because  of  the  required  control  of  the  stimulus  data.  Such 
control  cannot  be  achieved  through  the  use  of  the  acluci  system.  The  role  of 
the  real  world  system  in  the  laboratory  concept  was  played  by  a  computer  simula- 
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tion  of  a  soles  distributor  company.  The  SIMCO  simulation  was  used  because  of 
its  known  charocteristics.  The  outputs  of  SIMCO,  considered  as  standard  for  com¬ 
parison,  were  matched  with  the  outputs  of  the  DO^JA  model.  In  this  manner 
a  judgment  could  be  made  on  the  performance  of  the  DONA  model  and  the  method 
used  to  produce  the  model. 

The  Investigation 

The  investigation  was  centered  oround  the  feasibility  of  developing  the 
approach  and  converting  it  into  a  suitable  computer  program.  This  portion  of 
the  project  was  indeed  a  huge  undertaking.  The  DONA  approach  requires 
extensive  matrix  manipulations.  The  matrix  manipulation  would  not  be  such  a 
problem  except  for  two  very  peculiar  problems  ossociated  with  modeling  a  large 
man-organized  system.  To  appropriately  chorocterize  a  large  system,  many 
variables  need  be  considered.  An  example  of  this  is  the  relatively  simple  SIMCO 
simulation  which  contains  about  50  variables,  if  all  these  variables  were  con¬ 
sidered,  the  DONA  method  dictates  t'le  mathematical  operations  of  matrices 
with  a  size  100  by  100.  The  second  problem  that  arises  wos  one  that  was  com¬ 
pletely  unforeseei  in  the  early  stages  of  development.  It  appears  that  in 
quantitatively  characterizing  a  man-organized  system  many  redundancies  occur. 
This  results  in  linear  combinations  appearing  in  the  matrix  descriptions  of  the 
system.  The  DCNA  method  requires  inversion  of  the  previously  described  large 
matrices.  The  presence  of  linear  combinations  In  a  matrix  to  be  inverted  causes 
well  known  problems.  This  problem  was  resolved  by  removing  the  linear  combina- 
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Hons  and  fhereby  avoiding  matrix  singularities.  A  third  problem  which  arose 

was  due  to  computer  capacity  limitations.  The  system  description  matrices  con- 

9  -6 

tain  numbers  ranging  from  10  to  10  .  The  problem  was  solved  by  using  o 

"normalized"  form  of  the  matrix.  The  motrix  was  converted  to  intercorrelations 
which  range  from  +1  to  -1.  With  the  above  major  problems  resolved  the  computer 
program  was  complete. 

The  investigation  then  proceeded  into  the  validation  phase.  The  input 
or  stimulus  data  was  conceived  of  consisting  of  one  or  a  combination  of  impulses, 
step  functions,  si  nusoids,  and  random  noise.  A  comparis^'n  was  made  of  the 
outputs  of  SIMCO  and  DONA  model  when  the  same  stimulus  was  applied  to 
each  model . 


Results 

The  primary  objective  of  developing  the  DONA  approach  and  con  verting 
it  to  a  computer  program  wos  achieved.  This  program  is  shown  in  the  flow  chart 
in  Figure  4.2.  The  computer  listing  for  the  program  is  in  Appendix  D. 

The  validation  results  are  shown  in  the  computer  printouts  in  Appendix  E. 
It  will  be  noted  that  the  DONA  and  SIMCO  outputs  do  not  closely  compare. 

The  reason  for  this  is  explained  in  an  earlier  portion  of  the  same  printout.  The 
program  contains  a  unique  feature  of  evaluating  the  "framing"  of  the  problem.  It 
is  important  in  the  initial  stages  of  modeling  to  carefully  select  the  real  world 
variables  for  analysis.  It  is  highly  unlikely  thot  all  variables  may  be  entered 
into  any  model  of  a  man-organized  system.  Therefore,  a  careful  screening  action 
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most  be  performed  to  select  variables  thought  to  be  most  influential  in  the  per¬ 
formance  of  the  system.  The  screening  action  that  took  place  in  the  SIMCO 
liimulation  is  a  good  illustration  of  the  "froming"  requirements.  As  stated  earlier 
5IMCO  contained  about  50  variables.  After  analysis  of  the  significance  of 
each  variable,  output,  2  input,  and  5  state  variables  were  identified  for  use 
in  the  modeling  effort.  Then  the  computer  program  evaluates  the  contribution 
of  these  variables  toward  explairiing  the  variability  of  the  model  outputs,  This 
accounting  for  the  variability  in  the  dependent  variables  (outputs)  as  explained 
by  the  voriability  of  the  independent  variables  ^inputs)  is  reflected  in  the  co¬ 
efficients  of  correlation  in  the  multiple  regression  analysis.  The  square  of 
the  coefficients  may  be  considered  os  percentages  with  100%  meaning  perfect 
correlation  between  dependent  and  independent  variobles  and  0%  indicating  no 
correlation.  The  computer  results  show  that  the  correlations  range  from  0.4% 
to  58%.  Then  the  poor  comparison  of  model  outputs  is  to  be  expected.  But  whet 
this  Is  saying  is  that  there  have  not  been  either  enough  voriables  entered  in  the 
model  to  characterize  the  system  performance  or  the  ones  currently  under  con¬ 
sideration  ore  not  contributing  enough  to  explain  the  system  performance.  An 
illustratiori  of  this  is  an  experiment  actually  performed  to  try  to  increase  the 
"accounting  for  variability."  Prior  to  accomplishing  the  58%  that  variable 
had  given  results  of  about  23%  when  anothrrr  choice  of  quontities  had  been 
entered  in  the  model.  The  35%  increase  was  due  to  dropping  out  one  variable 
measured  from  SIMCO  and  substituting  another  one.  It  is  felt  that  further  work 
on  this  problem  would  yield  better  comporobil ity. 
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The  concept  of  modeling  through  the  state-space  approach  and  the  DONA 
methodology  looks  extremely  promising.  This  thesis  was  o  beginning  in  the 
development  of  the  concept  and  implementation  through  a  computer  program. 

It  is  recognized  that  many  areos  have  not  been  considered.  It  would  be  desir¬ 
able  to  explore  the  basic  characteristics  of  the  DONA  model  and  the  method 
used  to  produce  the  model.  Detailed  experimentation  in  the  laboratory  system 
would  provide  some  of  these  answers.  The  laboratory  system  was  >.onceived  to 
permit  control  of  the  variables  to  allow  discovery  of  these  basic  characteristics. 

A  systematic  application  of  the  above  alluded  to  types  of  stimulus  date  would 
lend  understanding  to  the  responsiveness  and  frequency  characteristics  of  the 
DONA  model. 

Further  effort  can  be  devoted  to  enriching  the  conputer  program  to  auto¬ 
matically  screen  the  system  data  to  detect  the  possibility  of  linear  combinations. 
Since  this  condition  causes  an  automatic  "no  solution"  in  the  DONA  methodology, 
this  area  should  get  priority  consideration.  It  is  believed  that  this  problem  can 
be  easily  overcome. 

Also  fr  r  further  study  is  the  area  of  comparison  in  the  laboratory  system. 
The  comparison  at  present  is  based  on  a  judgment  of  the  comporabil ity  of  two 
columns  of  figures.  A  statistical  means  should  be  odopted  to  permit  a  mere 
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This  statistical  procedure  lends  itself  to  an  analysis  involving  multiple  variables. 
Further  information  on  this  proposal  is  contained  in  a  paper  by  B.  B.  Lukens 
and  R.  A.  Brown  entitled  "Experimental  Design  Considerations  in  Validating 
a  Mv'thod  of  Modeling  a  Man-Organized  System, 

Explorotion  of  these  areas  would  add  tremendous  value  to  the  DONA 
concept  of  modeling  a  large  ma:; -organized  system.  A  significant  achievement 
would  be  realized  by  applying  this  concept  of  modeling  to  the  real  world  system. 


L. 


APPENDIX  A 

COMPUTER  PROGRAM  FOR  SIMCO 
AND  SIMCO  OUTPUT 
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w  fOH  SIMCO 
(.  PUUOMAH  StMCO 

OIMLNjION  AH(bl fPlPLlN(9l»00P<ltl»R(22)>A(22) 
mHtNblOU  S1LitV<22l 
l.OWWON/RA'M3/lHD«N5T*BT«l»NDlilRNiKiBPl 
r4STAMrsi2J>tba 
Call  inji 

1  F0HMAMI<4M  dttKLY  0£MAN0.F27.HiF2b.4/23H  BACKORDER  QUEUE  LENGTH, FI 

iB.i«,F^b.4/2bM  Inventory  at  eno  Of  MCeK,Fi6>4>f25,4/i7H  supply  back 
20RutK.F24.4.F25.4/20H  ACCOUNTS  RECElyABLEif2I.4»F2S.4/l3H  CASH  uN 
3HAND,F2e.<«>F2b.4/bH  SALES, f3b*4,F2S.4/24H  6R0SS  INCOME  FROM  SALES* 
4K1 /.«4.F2b.4/2UH  meekly  GROSS  PROFIT, F2I.it,F2S,4/18H  MEEKLY  NET  PRO 
Sf 1  I .►rf3.4.F^a.4/20H  UNITS  OF  LOST  SAL£S*F21 .4 ,F2S.4 I 

2  FoMRAl t//33A,eH  AVERAGE tlOXtlRH  STANDARD  DEVIATION) 

J  FCMMAl  l/’//JX,4HMEEK,bH  AUUM,bX,3HCR0,SK,4HAIEW,S)(,3HASB,7X  ,  3HTAh, 
14X,3HLUH,bX,3HTUS,«iX,3HGIS,7A,3HMGP,7X,3HWNP,4X,3HULS*SX,2HFD,3x,3 
2HPMU/ I 
MKITl  (h,3l 
(.  luniALUt  VALUES 

V.  UiITIAlWE  data  accumulation  values  at  2ER0 
Uo  4  ksi,22 

4  H(K*  :  U.O 

L  lulTlALUL  PIPELINE  FOR  THE  NEXI  0  MEEKS 
DO  S  M:),6 
S  PlPLlNlMtsbU.O 
PIPL1h( 71:0.0 
PJPLlN<8l:0.0 
Tpu  :  SO.O 
SU  :  15. 

UOSbu.O 
ULS  :  0.0 
AlEkSibU.U 
C(<g  :  0.0 
00  t  J  =  1*10 
t>  UUP(J):bO.O 
Fu:bU,U 
PRO  =  0.0 
SSF  :  O.U 

0  JiiITIAlUl  accounts  receivable  Fur  the  next  4  MEEKS 

DO  7  k:1«4 
7  ARIK):2S0UrU 

tamsiuoou.o 

Coh:400U.0 
L  HUN  FOR  luo  hEEkS 

UO  2b  JmKS=1,100 
CALL  HANUMIA.UI 

L  call  supply  10  uetermine  leap  time  for  units  orderep 
call  svjpply IPIPLIN«SR«U0,H1 
L  *  ill  CUSIuMtK  dalhOROERS  if  Any 
AIUM  :  AlLM  *  SR 
U  IChO  -  AlbMl  H>A,<) 
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8  us  :  CUU 

AAIUW  =  AIB«  •  CBO 
Cbu  =  O.u 

ULS  =  O.U 

&0  TU  10 

V  US  =  AlBW 
AA1B«  =  0.0 

ULS  =  (CUd  -  AlBMl  •  .20 

CUU  =  <CUU  >  AlBUl  •  .80 

lU  N  ::  CUU 
CUU  :  N 
N  =  ULS 
ULS  =  N 

C  call  UtMANU  TO  ULTERMlNi  OCMANO  HOH  THIS  WCCK 
CALL  UtM AND  I AUOW • PMO . TMU . SO • SSF  »  A ) 

L  uLTCHMlut  TOTAL  SALES  DURING  THE  MEEK 
IT  (AuUH  -  AAIBH)  11«11.12 
11  TUS  =  US  «  AUUW 
GO  TO  13 

1^  TUS  :  US  *  AAIBW 
«.  ULTEHNlNL  BACKOHUER  QUEUE  LENGTH 
CUU  :  CriO  «  (AOOh  -  AAIBW I 

C  uLTEHNllit  GROSS  INCOME.GROSS  PROTIT.  AND  ENDING  INVENTORY 
13  GtS  =  lOU.  •  (1.0  -  PMOt  «  TUS 
WbP  =  61S  -  SO.O  •  TUS 
A|L«  :  AAlBw  •  AOOw 
IF  (AlLh)  U.IS.IS 
lu  Altw  =  O.U 

IS  continue 
lb  FCS2UUQ. 

I.  UlTERRlNt  CASH  ON  HANU 

17  COH  s  COH  4  AH(U  4  .S  *  GIS  -  SM  •  SO.  -  FC 
1-1.2/32.)  •  SU.  *  AIBW  -  (.0tt/S2.)  •  TAR 
C  UCUATL  ACCOUNTS  HECElVABLE  AND  DETERMINE  TOTAL  A.H. 

DC  lU  J  =  1.3 
lb  AH(J)  =  AH(u41) 

AN(«|)  =  .S  •  GIS 

TAM  S  ARID  *  AR(2)  *  AR(3)  *  AR(b) 
u  uLTEHMlNt  weekly  net  PROFIT 

WNP  :  WGH  -  (.2/S2.l«S0.  •  AIBW  -  (.08/52.)  •  TAR  -  EC 
IF  (AUOW  -.7U  •  FO)  19.21.21 

C  FUHECAST  UtMANO  FOR  THE  COMING  WlEK  AND  DETERMINE  IF  A  PRICE  MARKDUmN/ 
C  IS  NECESSARY 

19  CALL  FORCAS(UOP.FO.ADDW> 

IF  (AIEW  4  PlPLINdl  -  1.5  •  FDI  22.22.20 

20  (>Mu  :  O.IS 
SSF  s  U.US 
GO  TU  23 

21  CALL  FORCAS  (OOP. FO. AOOW) 

22  PMU  :  U.O 
SSF  s  0.0 

23  CONTINUE 

«.  call  ORUEh  to  ULTERMInE  the  number  OF  UNITS  TO  BE  ORDERED 
CALL  UNDER  (PIPLlN.UO.TuS.AlLW.FO.ASB.DI.nsB) 

WRITE  (b.2<«)  JWKS.AODW.CUO.AXEW.ASB.TAR.COH.TUS.GIS.MGP.WNP.ULS.FD 
•  .PMD 

2<4  FORMAT (I7.Fb.U.Fa.2.2F9,2.F11.2.F12.2.Fb.fl.3F10.2.F6.1.F6.0.FS.2) 

V.  ACCUMULATE  SIGNIFICANT  DATA 
R ( 1 )  =  N ( I )  4  «D0w 


H(^)  S  H(;<  *  A00W««2 
S  Hljl  *  CHU 
Ntm  S  HU)  *  CBQ«*2 
HCjI  =  H<St  ♦  *IC«( 

Hl0>  •  H(6t  *  AICM««2 
H(7)  s  MIT*  ♦  MSB 
RIB*  S  RIB*  *  *Sa««2 
R(9)  S  RI9>  ♦  TAR 
HIlO)  =  Milo*  ♦  IAH*«2 

HUH  =  Hill*  ♦  COH 

HU2*  =  HIU*  ♦  C0M**2 

HU3*  :  MllA*  *  TUS 

RU«>  :  HUM*  «  TUS«*2 

HU&*  =  M(lb)  *  Olb 

Hilo*  =  Mllb*  *  aiS4«2 

hut*  s  mIITI  ♦  rBR 

RUB*  =  HUB*  ♦  Ba»»**2 

MU4I  s  MIlV*  ♦  HNf 

HIBOl  s  HI20*  ♦  MNH««2 

HIBI*  =  MIBl)  *  Utb 

|f^»U*^l‘j"l.tSS*^  ImAn'^100**RCTURN  to  T*C  lOe  OF  THC  LOOP 

re  sieiriCNT  o.,. 

UO  Ba  UtrBlrB 

siotwn *:‘»CiNi I iRi I *1  >»Hn I ••2/ioi>»>/m*> 

I  MlltT  COLUMII  HLAL1N6S 
•HlIL  IB<2) 

I-hJI.T  Ht'jULTS 

•MlTt  IBrll  IAlI*#!»TOEVUI.Ul#2l*2* 

STOP 

ii.0 

f  UH  01  M*l.u  .  .  ,,, 

SOOHOoTlt.t  oLMANUIAOO*l*PMO«TMil»SU»SSF*V* 
c  S  IMO  »  So  •  V 
1*  10*  UBr^ 

1  0  :  UrO 

2  1*  IPMO  -  Oris*  IrHrH 
^  SSP  s  O.ll 

(iO  TO  5 
0  SSP  :  O.us 
.  AOOm  :  0  *  UrO  *  SSP* 

N  S  AoOh 
AbU«  :  N 
HE  Tun IV 
LMO 

SOBHOUTINE  SUPPLY  |PlPHNtSR»OOt A* 

UlMtNSXON  P1PLINI9I 
OtTEHMjNL  SHIPMENTS  RECEIVED 
SN  S  PIPLINU) 

UPDATE  SUPPLY 
00  2  ASl.T 

B  PlPLlNlA*  s  PIPLINIK*!* 

PlPLlHlB)  s  U.O 
oETEHMliJL  LEAD  TIRE 
MsiUuU. 
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J  t  .T  s  <♦ 

GO  TO  11 
4  IT  IM*24I 
b  UT  =  t> 

GO  TO  11 

t,  If  jM-7m  7«7<8 
7  (LT  s  u 
GO  TO  11 

d  IT  CM-SV)  9*V»10 


GO  TO  11 

.  tr«:l;;«L°.l"WO»Kl.T|-  .«.s  SU«-L»  WOtKO  «  «.US  WOtMO  THIS  «l,» 

PieCl.>.ILtl  s  HIPLINILTI  ♦  00 
«tIUHn 
CNU 

*  KOM  FUMCA^j  _  , 

bUMHOUTINC  TO><C*S<00P*P0**00»* 

Uli^CNblON  UUP(11>  ,  „ 

C  uOGNOATt  OEMANU  FOU  TM€  UAST  TEN  liEfcHS 
00  1  lbsl*9 
L  =  11  -  10 
I  UOP*l-t  =  UOPit-l> 

OOP*  It  =  *OUli 

(.  CALCULAIt  FOHEtAbTEO 

FUSIOOP*  ll^toPt2»^OOP<SH/a«0 

HE  ruHG 

LNU 

*  SouTint  GKUE«IPIPL1N.OO.TOS.A1E«.FO»ASB.OI.DSB» 

OlMtNblOr4  PlPUINim 
1.  CALCULArt.  OtSlUtU  1N0ENT0R7 
Ol  S  J.O  •  TO 

C  (.ACCOUATt  SOPfLt  bACAOPOEM 
«SU  i  U.U 
00  1  H  =  1»» 

I  A-J  S  ASO  ♦  P1PLJN(N» 
t  oAtCULATt  OEblHtU  SUPPLY  eACROPDLB 

So^/'tS;"**.^^  «0I  -  *IEG  ♦  OSB  -  ASO) 
n:uo 

UOSN 
If  (out 
i  00  s  U<(t 
3  CGhTIhuE 

Ml  TOPI* 

ttiU 

.  » Ok  KANUM.hanum 

bOOHOuTlNE  KANOPtOtViOT *HN02I 
tl.AM0.4/HANO/JHO«NbTA«T.MN01*i»N.NOPS 
IKN- l«bS< |Hu**2*NOPS**2t 

iKo=<iHN-aKN/iooooooooo»*io«oMooooo»/ioon 

If  » IHo.tO.O) 1H0=1 

NOPSSHOPb»l 

HNU1:kNU2 

MLTUPK 

ENU 

,  FoK  irdTtlMT 

SOBhOoTlrJl  INIT  „  ^  _ 

C0PP0i4/MAh0/ IPOfNSTAPT  fBNOl ♦ *MN*NOPS 

IMUSNSTAKT 

MNOlsFLOATtNbTAPT) *0.000001 

r40PSS4743b 

RETUHU 

LNU 

•  «ul  SIMCO 
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f»r»rncrt'OCo4‘i 
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APPENDIX  B 
MODIFIED  SIMCO 
COMPUTER  USTING 
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b  FOH  SIMCO 
C  PRO&MAM  SIMCO 

dimension  AM(S)>ARl(SI*PlPt.lN(«>*f>lPLlE(9>iOOP(U)>tjO«>i<ll))V<i01) 
i>viiiui)>xM<ioo)tXMiaoo) 

COMMON/R«NO/lHOfNSTART*RNOt«IRN*NOPS 

NST«RTsX23‘*S6 

CALC  INIT 

3  FORMAT  {///mi itHNECK •  HHAOOM •  HX •  5HA00H 1  ( 6X •  3HSAt >6X » 3HF 01  •  7X •  4HT US  1 
*•  6Xi3HC0HiSXi3HTUS)«Xt3H61SiTX*3HW6P«7Xr3HHNP«4Xf3HULS>5X»2HFD«3X 
••2HZT>3X*lH2i3X»2HMT>4X»3H  SM/I 
taHlTL(6»3> 

(.  IMTIALUL  VALULS 

iriiTlALUL  PIPELINE  FOR  THE  NEXT  R  MEEKS 
00  S  H=I.6 
S  PIPLINIMISSU.O 
00  SO 

SU  PlPLiL(M):SU<U 
PIPL1NI7I:0.U 
P|PLlt(7ISO.U 
PlPLlNiaisO.U 
PlPLlt(8i:U.U 
TMU  =  SU>U 
TN<U1=SU.U 
S0=2U.U 
SCUld.U 
uossu.u 
uoi:su>o 

ULS  =  U.U 
ul:31  =  u<0 
AILmSISU.U 
AlLoUlSU.O 
Ci<g  =  u.o 
cooiso.u 
uo  b  o  :  1 • 1 U 
o  UCPtJtsSU.O 
00  bU  j:l<lu 
uU  U(iPl  *J>:bO.U 
FOSSU.U 
F01SSu>ll 
Pt'U  5  0.0 
P^UlSU  .0 
SSF  S  U.U 
SSFISU.O 

c  IMTIAlUL  accounts  RECCIVARlE  Fur  the  next  4  WEEKS 
00  7  b-l>4 
7  AK\ms2bUU.U 
OC  7U  K=i.4 
/U  AKKis  tsZbOO.U 
TAhsIuUOU.O 
TAMlSlUOOO.u 
COH=WuUO.U 
COHt=4000.0 

<.  lutllAcUL  PERSONNEL  ACTIONS 

SUMM  =  0.0 
SUMMISO.O 
T(l)S3.0 
T1 ( I  1=3.0 
XM  1  1=0.0 

XMiiitro.o 
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C  RUN  fOR  100  MEEkS 

00  2:>  JHKssi.loo 
CALL  HANOMUtOI 

C  CALL  SUPPLY  TO  UETERMINE  LEAO  TIME  FOR  UNITS  ORDERED 
CALL  SUPPLTiPlPLlNfSRiUQiBl 

C  FILL  CUSTOMER  BACKORDERS  IF  ANY 
AIBH  =  AIEm  ♦  SR 
IF  <COO  -  AlbH)  6ta*9 
0  US  :  CUO 

AAlBa  s  AIBk  -  CUO 
CBU  s  0.0 

ULS  s  0.0 

60  TO  10 

0  US  s  AIUh 
AAlUa  :  0.0 

ULS  -  ICliO  -  AIRmI  •  .20 

CbU  =  (CUO  -  AIRhI  •  .00 

lu  U  Z  Lou 
CUU  =  U 
N  S  UlS 
ULS  =  N 

L  CALL  ULMAuU  TO  UtlERMiNfc  DEMAND  FOR  THIS  WEEK 
CALL  uLMANU(AUOh.PMO«TMu>SO«^SF*A) 

L  OLTLHMlf4t  total  SALES  UUMIN6  THE  MEEK 
IF  (AuU«  -  AAlHMt  ll«Ufl2 

11  TUS  =  US  «  AUUW 
GO  TO  n 

12  TliS  s  HS  «  aAIUM 

L  uLTEMMIiL  UACHOHoER  UUEUC  LENGTH 
CUU  :  CHG  *  (AOOm  -  AAlRw) 

L  ULTLHA’lnL  GROSS  INCOME.GROSS  PROFIT.  AND  ENDING  INVENTORY 
11  GIS  z  IQU.  •  (1.0  -  PMO)  •  ToS 
nGP  z  otS  -  SU.O  •  TUS 
A1L«  =  AAIHm  -  AUUm 
IF  (AlEN)  U.IS.IS 
1<*  Albs  S  0.0 
IS  CONTINUE 
lb  FCS77/.0 

L  UETENMlNb  CASH  01.  HAND 

17  COHSCUHFAM (114 .S*G1S-SN«S0.>( .2/S2. US0.«AIHH-( .0A/S2. I •TAR 

L  UPUATE  ACCOUNTS  NLCCIvABLE  AND  OlTERMINE  TOTAL  A.R. 

UU  lb  J  s  1.1 
lb  AF<(JI  Z  AH(o41l 
AN(iil  S  .S  •  GlS 

TAM  S  AHdl  4  AR(2I  ♦  AR(1I  4  ARIbl 

L  ULTLMMINL  mELALT  net  PROFIT 

•(iP  s  «6P  -  (.2/b2.)*S0.  «  AlHu  >  l.se/S2.>  •  TAR 
IF  (AuOm  *.70  •  FOI  19.21.21 

C  FORECAST  UEMANU  FOR  THE  COMING  wt.E((  AND  DETEHMINC  IF  A  PRICE  markuo.i./ 

C  IS  NECESSARY 

IV  CALL  FOMCAS(UUP.FU>AOOMI 

IF  (AiEH  4  F’1PL1N(1>  -  l.S  *  FD>  22.22.20 

20  PMU  z  U.IS 
SSF  s  O.OS 
GO  TO  21 

21  call  FORCAS  (UDP.Fd.AODWI 

22  PIAO  s  0.0 
SSF  s  0.0 

21  continue 


C  call  OHOfcK  TO  OtTEMMlNC  THE  NUMBiR  OF  UNITS  TO  BE  OftOEREO 
CALL  ORDER  (PlPLlN*UO«TUS> AItM*FO» ASB*OI >DSB) 

C  LNO  Of  PNUOUCT  !•  8C6IN  PRODUCT  2 
call  HAN0M(A«B» 

CALL  t>UPPLl<PIPLIE*SRl>UOIiBI 
AIHMlsAieiil«SHl 
If  «CULl«A|8kill  80>60>«0 
HO  BSlsCUOl 

AAlU«l=AI(tMl>CB01 
CHUUU.O 
UL!>i:u*U 
00  TO  lOU 
00  BSI=A1BMI 
AAlUMlSU.U 

Ul^UICHul-AlUMl  )*>20 
ChOi:<CHUl-AlHHl>«.60 
lUO  NSCUOl 
CbUi:i< 

NSULSI 

ULblsi* 

call  uLmnu  ( AU0«l>PM01>TMUl«SUl«SSFlf A) 

If  (AWOltl-AAlHIIll  Il0«llb>t2u 
iiu  rusisusi*AOu«i 
60  10  ISO 

1^0  TU!>13HSl*AAlt>Hl 

CbOULBOl*  ( AUUWI-AAIBWI I 
tJO  61SI:IOO.*IUU«PM01I*TUS1 
«bPlsolSl-Sb.U«TUSI 
AlLkl=AAlU«iI>AOOMl 
IF  (AiLall  lH0»lSUilS0 
IbU  AILblsO.U 
lt>0  CONTINUE 
ItiO  fCl:7/7.u 

1 70  C0Mlsc0tU*AHl<n«.b*CISl-SRI*5U.>l. 2/52.  >*S0.»A1HH1>I. OR/52.  >*TA><1 
liu  IbO 

loU  Am1(J):AH1(.<*1I 
AhKbIs.bvGISl 

TANIsAHUI  )«AH1  (2l«AHll3l«ARiU) 
aKHi:«6Pl>( .2/52. 1*50. •AIHOl>(.0i>/52.) ATARI 
IF  iAuU«I«.70*f01l  190>2t0«2i0 
I'lO  call  f  ORLST  IOuPl>F01«AOO»flI 

IF  lAieMlAPIHLlEUI'I.BAFOIl  220«220t200 
^uo  PPoiso.lb 
bSf i:u.5 
OC  To  230 

210  call  f OMCbTIUUPI«FUl«AOOkII 
220  PMUiSu.O 
Sbf ISU.O 
230  COHTIHUL 

CALL  ON0NlPIPLie>U01>YUSI>AlLMl.FUl«ASBlfr)U»0SBII 

CALL  »'CNSI<IL  lfp.FUl.V«VI.MT.2«2l»2TrSALtSALl*SUMM»SUMMl>JWKS»M.Nl) 
L  Ut'OATL  COMiuL&.mNP  WITH  PERSONNEc  CONSIDERATIONS 
COHStoH-FC-bAL 
COHISCOHI-FCI-SALI 
WMf>T=  ( WNH  AWUPi » -f  C-f  C  1-SAL-SaL  I 
J»KS2=JwKbAl 

af<ITL  <bi241  JHKb2<AU0w|.F0I.C0H»6IS*NNPT.FU 
/<•  fOHMAr </17.tiA.FB.2.9X.F9.2.IlX.F12.2»6X.f 10.2.lnX.FI0.2.6X.Fh.U) 
ahlTL  (6.35>  A00«>bALfTUSI.TuS>W6P»ULSt2T.2»MT.SH 
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Jb  FOHMAI  <7X»l-b.0»«X«F9.2*9)(tF11.2*12X«r6.0#10X*fl0.2»10X>F6.1»b«»r 
*!>.2.2XtF3.1>2X«Xa»2X«FS.2l 
CONTINUE 
STOP 
ENJ 

»  FOH  OCMANu 

SUUKUUTINE  UlMANO(AOOW*PMOfTMU»SO«SSF*V> 

L)  :  IMU  4  SU  «  V 
IF  <0) 

1  U  =  U.U 

<!  IF  IHmU  -  U.lbl  3»4>A 

3  SSF  =  U.U 
60  TU  b 

4  SSF  :  U.Ub 

b  ALU*  =  0  •  ll.n  4  SSF I 
N  =  At>UM 
AUUk  =  N 
HLTuki. 

LFaU 

..  FUM  SOPPLf 

Sl'UHOoTINk  SUt'PLr  (PIPLIN.SPtllOtXF 
OlMClijlOh  F>1PL1N(9> 

S4  :  PlPLlNdt 
'JO  2  4:1  •  7 

<  F'lF'Ll  .14  I  =  PiPLlN(K4il 
PlF<Ll>.lMi  :  U.O 
MiA«lUO. 

IF  IM-91  ititH 
i  U  Z  >* 
bU  1U  11 

4  IF  (M-24i  S«b46 
b  LT  S  b 
00  lu  11 

«i  IF  <M-7m  7./*A 
7  LT  s  « 

60  TO  11 

b  IF  (M-a9)  9f9*10 
9  OT  s  7 
60  TO  11 
lU  LT  S  M 

11  PlPLll4<L1l  s  HIPLlNtLTI  4  UO 
NE  rUMl4 
Er.u 

«  FjH  FONCAs 

SUrtHUuririL  FUMCAS<UOP>FU«AOOat 
uiMiNblor^  uoPiil) 

UU  1  iu:l«9 
L  S  11  -  IB 
1  UOPItl  s  U0H(L-1I 
UOPtlt  r  AOU« 

F0=IUuF>ll)4bUP(2i4U0PI3>  )/3.U 

HtTuMi. 

bt.U 

*  FOO  OHUI.m 

SUbOOijTlHt  OMUEmPlPLIFI.UOfTUSfAlEW.FO.ASB.DlfOSbl 
UlMCNblOn  P1PLINI9) 

01  :  3.U  •  FU 
ASU  z  U.U 
UC  1  M  Z  1 9  D 


1  ASH  =  ASH  *  PlPLlN(N) 
usu  :  6.U  •  FU 

UO  :  tUS  *  .b  «  <01  >  AlCii  *  OSB  •>  *SB) 

N=UO 

UO=N 

l<  ( Uu )  / •  3> ^ 
t:  UO  :  u>U 
b  Cui«1  IWUL 
wL  ruHij 
tNU 

..1  I  FOH  HANliPiHANUM 

SUUMUoTlNt  NAnUM(ULWIOT  it4N02> 

CCMMUii/HAr40/  IHOiNSTAHTiRNOl  •  IPNtNOPS 
1KN=1AUSI IMu*«2«NOPS**2) 

ihu:i iHN-i tHN/iooouooooQ)«iounoooooo)/iooo 

U ( iKu.tU.O) IHOSl 

NUPS=i«OPS*l 

HF4U2:<  LUAT  (  IMU I  *0.000001 

Ut  V 1 0  T :SOH  T  < -2 . • AlUO  <  HNO  U ) aCOS  t  b . 2BS 1 BS4*HN02 ) 

HNUtSHNU2 

Hk  rUHU 

tt.U 

*11  FOP  Ir.IlitUlT 

SuUHUgTiNi.  INtT 

(.OMMUx/RANU/  lHOtN<_TART  (PNUl  •  iRNrNbPS 
IPJSN^TAHT 

PliUl=<-'l.OAT  (NSI  APT)  *0.000001 

N0PS=**7P<>b 

HI  rUHN 

tNU 

FUH  ULMNU 

SOBHUoTINL  oLhND  ( AUOMl tPHOl t TMUl fSOI <SSP t i V > 

0  s  TMU1*S01«V 

U  <U)  1>2>2 

1  USJ.U 

2  If  <PMU1-U.1SI 

1  SSFlsu.O 
OO  TU  S 

4  SSFlsu.Ob 

5  AUUi»UO*<  l.H*SSPlt 
NSAUU*! 

A0U4UN 

HITUHh 

LNU 

.  KJH  SUHPLl 

SUHHwoTluL  SUPPLt  <PlPLlL>SHi.U01tXI 
UlMLnblON  PIPLIE  (91 
SHlcPlPLUdl 
UO  2  a:1i7 

2  PIPLU  IAI=PlPLie  (K*l) 

PIHLlt.lai=0.0 

M=X*1U0. 

If  <P-9»  i.J*4 
i  LT  5  4 

00  TU  11 

4  ]f  (M-24)  S.bib 

t>  LT  =  b 
00  To  11 

b  IF  IM-741  7.  he 
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7  LT  =  6 
(iU  TO  11 

a  IF  (M-a9)  9«9>10 

9  LT  =  V 
bO  TU  11 

10  LT=«  . 

11  PIPLIi.lLTiSF'iF’LI£<LTl*U01 
HCTUHU 

ENU 

*,  FOH  F0»C!»1  , 

SUBHOoTINE  FOKCST  (DOPl *POl • AOOBU 
UlMENblON  OLlPl<ll> 

00  1  1H=1>9 
L:11*1U 

1  OOPl  (DsuOPl  St-l  I 

OOPUl»s*DD»l  „ 

F01slt)UPUl)^00PH*l*0DPl<5)  »/3»0 


METUHl’l 

ENU 

**  SUBPOuTINt  OKUR  <PIPUlE.U0l»TU&l»AIEl(l*fOl**SB1.0Il*OS8l) 
OlMENSlUN  PlHtlE  <91 
U11:J<U*FU1 
ASblsUtO 
UO  1  NSl'A 

1  A<,U1SA&B1«P1PL1ECN) 

USB1So.0*F01 

U01  =  TOil*.5*lUU-AlEBl«0SUl-ASBn 

NsUOl 


UOiSN 
IF  (UOlt 

2  UOISU.O 

3  CUNTll<Ut 
Ht  tUHN 
ENU 

mu  PtKSNl. 

SUBRUUTINE  PLKSNL 


(F0»F01*T*Tl»MT*2»Zl»ZT*SAL»SALl*SUHM*5UMMl.d*AS 


OlBiNblON  T*  lull  •  n  U01».  XMUOOl  tXMi  1100  J»XMT<  100) 
T<J*(AS*1)sFO/20. 

TH  Ji»Kb*l)SF01/20. 

L2=V«JBKi*l 1*0 *9999 
tlsTt J«KSI«U*9999 
MSkZ'Ll 

L12STI<U»IKS*1»  *0.9999 
Lll  s  yi< JBK>)*0.9999 
Ml=tl<*tll 
X»)(J«aS*1>3M 
XMl<JaKS*l)SMl 
SUMM=bUMM*  XN ( JMKS  > 

SUMM1SSUMM1*XB1(JNKS) 

2SSU0IM«V  1 1 ) 

21=SUHM1*V1(1) 

MTSMl+M 

*tb  XMTIU«KS)=MT 
ZT521*Z 

IF  (XM< JaKS) )1*3*9 
1  SAt.r2*192.-XM(JWKS)«2t«192. 

GO  TO  A 


i 

4 

b 

b 

7 

XUT 


b«Ls/*l92.  .  ,  ^ 

IF  UHUJbKS))  bt6i6 

S*Ul=2X*l«.-*Xl « JbUS)**.  < 

GO  7U  7 

SALl=2l«X92i 

RC.TUMN 

END 

&1MCO 


Kll'i 

FIN 


-92  • 


APPENOIX  C 


MULTIVARIATE  UNEAR  REGRESSION  ANALYSIS 


77 


Multivariate  Linear  Regression  Analysis 


15 


The  problem  is  how  to  estimate  a  variable  or  set  of  variables  of  an 
unknown  universe  given  o  sample  set  of  data  from  that  universe.  The  method 
to  be  used  will  be  the  least  square  method  using  a  linear  model.  It  is  also 
assumed  that  the  errors  about  the  regression  are  normally  distributed.  This 
distribution  is  assumed  homoscedastic;  i.e.,  the  scatter  is  uniform  along 
the  line  of  regression. 

The  mathematical  model  for  the  universe  plane  of  regression  of 

cn  X-,  X.,, . .  X  is  represented  by  the  equation 
/  o  p 


X'  =a'  „  +b'  ,  _  X,+b'.,„  X.+  (1) 

1r  1,23. ..p  12.3. ..p  2  13.2. ..p  3 


. .  +  b'  ,  X 

lp.2o. , .p  -  1  p 


where  p  is  the  number  of  independent  voriables.  The  equation  for  the  esti- 
mote  of  the  universe  is 


X  =a,__  +b«„  X.,+b,_., 

Ir  i.23,..p  12.3. ..p  2  13.2. ..p  3 


(2) 


...+b,  „  ,X 

Ip. 23. . .p  -  1  p 


with  the  or ig if  at  X,  =  X^  =  X-,  -  •  •  -  =  X  =0. 

12  3  p 

The  primes  on  the  variubl-'s  indicote  the  universe  variables  and  are 
estimated  by  the  same  symbol  without  the  prime.  The  subscripts  are  divided  by 


a  dot  to  indicate  the  dependent  variable  to  the  left  and  the  independent 

variable  to  the  right  of  the  dot.  The  constant  a,  shows  that  the  depen- 

1 .23. . .p 

dent  variable  is  X  and  X.,,  X«.  ...X  ore  the  independent  variables.  In 
I  zap 

the  case  of  b^2  3  p  ^2  linked  together  for  the  dependent 

variables  and  X2...X^  are  the  independent  variables. 

The  method  of  least  squares  requires  rhe  coefficients  be  determined 
from  the  condition  that 


I(Xi-X,^  =MX,.  0,^23.. .p- 


*‘'12.3. ..p  ^2  ■  ...‘ 
Ip. 23. .  .p  -  1  p' 


be  a  minimum  for  all  sample  points. 

This  is  done  by  taking  the  first  derivative  and  setting  equal  to  zero. 

2 

Letting  (3)  equal  I  ,  then 

V  I  .z J . .  .p 


v^ 


1.23. ..f 


1.23...P 


1.23. ..f 


12.3. ..p 


*1p.2. . .p  -  1 


The  following  equations  are  the  result  of  applying  the  minimization 


ilihMUllftillfiillii 


T 


process. 
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no 


1.23. 


.p  ■^'^12.3. 


SX.,  + 
•P  2 


+  b 


Ip, 23. 


X  =2X, 
P  1 


'l.23...p^^2  ■^^12.3., 


2X/  +  ...  +b,  ,  SX-  X  =IX,  X_ 

.p  2  Ip. 23. .  .p  -  1  2  p  12 


1.23. 


IX  +b,- 
•P  P  >2. 


3.. 


IX.,  X  + 

.p  2  p 


+  b 


Ip. 23, 


(5) 

IX  ^  =IX, X 
P  >  P 


These  are  the  least-square  or  "normal"  equotions  for  estimating  the 
universe  regression.  Another  form  of  these  equations  can  be  obtained  by 
shifting  the  origin  from  0  to  the  mean  value  of  the  X 's  (X). 

The  resulting  equations  are 


-3  ^ +«-*+b,  1 

u,3...p  Z  1p.23...p  -  1  2  p  12 


(6) 


b./,..  Ix.,x  ^...^b_  —  Ix.x 

I2,i...p  2  p  1p.23...p-l  p  I  p 


where  the  origin  is  at  X, ,  X„, . . ,  ,X  . 

I  Z  p 

The  lower  case  x's  are  deviations  from  the  mean.  To  solve  equations 
(6)  for  the  "b"  coefficients  it  is  more  convenient  to  write  the  equations  in 
matrix  form  and  solve  by  use  of  linear  algebra.  The  following  matrices  are 


defined  as; 
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Then  the  matrix  equation  is 

AB  =K 


The  solution  to  equation  (10)  is 

B  =  a"V 


(8) 


(9) 


(10) 


(11) 


This  analysis  can  be  extended 


to  more  thon  one  dependent  variable. 
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The  same  procedure  may  be  used  wifh  other  dependent  variables  by  bringing 
in  new  "K"  matrices.  The  results  ore  the  new  regression  coefficients,  i.e. ,  new 
"B”  matrices.  A  compact  notation  con  result  if  the  "K"  matrix  as  defined  In 
equation  (10)  as  a  column  vector  is  augmented  with  more  coli^mn  vectors,  one  for 
each  dependent  variable.  Then  the  “K"  matrix  has  the  number  of  rows  equal 
to  the  number  of  independent  voilobles  and  the  number  of  columns  equal  to 
the  number  of  dependent  voriables.  The  "B"  matrix  has  the  number  of  rows 
equal  to  the  number  of  independent  variables  and  the  number  of  columns  equal 
to  the  number  of  dependent  variables.  The  "6"  matrix  will  have  the  same 
order  as  the  "K"  matrix.  The  coefficients  in  the  "B"  matrix  correspond  column 


by  column  to  the  dependent  variables  In  the  "K"  matrix. 

2 

The  standard  error  of  estimate  can  be  easily  determined  using  ^v  , 

I  •  Z'i  t  •  ej! 

which  is  equol  to  equation  (3).  The  unbiased  estimate  of  the  square  of  the 


universe  standard  error  estimate  is 


I .23. . .p  n  - p 


(12) 


where  n  is  the  amount  of  data  samples  and  p-lis  the  number  of  independent 
varables. 


The  squared  sample  multiple  correlation  coefficient  is 


The  squared  sample  multiple  correlation  coefficient  may  be  viewed  as  mea- 
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suring  tfie  relative  amount  of  variation  in  the  dependent  vorioble  that  is  "ex¬ 
plained"  by  the  independent  variable.  The  coefficient  of  correlotion, 

I  .£<}. .  .p, 

ranges  from  -1  to  +1  vrhere  a  -1  indicates  a  perfect  negative  correlation 
and  o  +1  indicates  a  perfect  positive  correlotion  between  the  dependent  and 
independent  variables.  Obviously,  a  zero  coefficient  of  correlation  says 
that  the  variables  ore  uncorrelated. 


APPENDIX  D 

COMPLETE  LABORATORY  SYSTEM  COMPUTER  USTING 
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DEHNITION  OF  VARIABLE  NAMES 
USED  IN  THE  LABORATORY  SYSTEM  COMPUTER  PROGRAM 


Main  Program 

The  main  program  of  the  laboratory  system  contains  three  major  subroutines 
”  SIMCO,  DONAFT,  and  DONA  —  as  shown  in  the  flow  diagram  in  Figure 
4.2.  There  are  three  arrays  of  variables  that  ore  common  block  to  each  of  the 
mojor  subroutines.  The  parameter  names  ore  also  defined. 


Varioble 
INDEP 
I  OB 

IPSI 

IXXX 

lYYY 

L 

MAX 


psi  a,j) 


X(1,J) 
YY  (1,J) 


Definition 


Number  of  independent  variables  In  DONA  model. 

Number  of  weeks  in  SIMCO  run,  number  of 
observations. 

Number  of  SIMCO  stote  variables. 

Number  of  SIMCO  inputs. 

Number  of  SIMCO  outputs. 

Same  as  INDEP. 

Sum  of  independent  ond  dependent  variables 
in  DONA  model.  (Number  of  independent 
variables  =  number  of  dependent  variables  in 
DONA  model). 

An  array  of  the  SIMCO  state  voriables  where 
the  "I"  index  is  the  time  variable  and  the 
"J"  Index  is  the  particular  dimension  of  the 
state  variables.  For  example,  PSI  (1,2)  is 
the  second  state  varioble  in  the  first  week. 

An  array  of  the  SIMCO  input  variables  with 
the  "1"  and  "J"  as  described  in  PSI  (1,J)  above. 

An  array  of  the  SIMCO  output  variables  with 
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the  "1"  and  "J”  as  described  in  PSI  (1,J) 
obove. 

SlMCO  Subroutine 

The  SlMCO  subroutine  and  its  internal  subroutines  are  shown  in  the 
flow  diagram  in  Figure  2.2.  The  internal  suboutines  to  the  SlMCO  subroutine 
Ore: 

DEMAND 

SUPPLY 

FORCAS 

ORDER 

RAN  DM 

INIT 

PERSNL 

The  following  is  a  list  of  the  variable  nomes  used  in  the  SlMCO  subroutine 
and  Its  internol  subroutines.  When  o  "1"  is  used  os  a  suffix  on  the  variable 
name,  the  variable  name  applies  to  product  two  corresponding  to  the  same 
variable  name  as  for  product  one;  eg.,  ADDWl  is  the  actual  demand  during 
the  week  for  product  two  and  ADDW  is  the  octual  demond  during  the  week 
for  product  one.  The  letter  "S"  is  used  as  a  suffix  to  o  varioble  name  to  change 
the  original  variable  name  to  a  new  variable  with  subscripts;  eg. ,  ADDW  to 
ADDWS(]).  * 

Variable  Definition 


A 


See  DEVIOT 


Variable 


Definition 


-AAIBW 

ADDW 

AIBW 

AIEW 
AR  (1) 
ASB 
B 

BS 

CBQ 

COH 

CS 

DDP  (I) 
DEVIOT 

DI 

DSB 

FC 

FD 

GIS 

M 

MT 


Adjusted  actual  inventory  beginning  new 
week. 

Actual  demand  during  the  week. 

Actual  inventory  at  the  beginning  of  the 
week . 

Actual  inventory  at  the  end  of  the  week. 
Accounts  receivoble  in  dollars. 

Actual  supply  backorder. 

See  RND2. 

Backorder  sales  at  the  beginning  of  the  week. 
Customer  backorder  queue. 

Cash  on  hand. 

Cash  sales. 

Demand  during  the  past  I  week  . 

Standard  normal  distributed  number  (mean  =  0 
standard  deviation  =  1). 

Desired  inventory. 

Desired  supply  backorder. 

Fixed  costs  in  dollars. 

Forecasted  demand. 

Gross  income  from  the  week’s  soles. 

Change  in  number  of  salesmen  for  coming 
week . 


Total  change  in  number  of  salesmen,  product 
1  and  2. 
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Variable 
N  START 

PIPUE  (K) 

PI  PUN  (K) 

PMD 

RND2 

SAL 

SD 

SR 

SSF 

SUMM 

TAR  (1) 

TMU 

TUS 

ULS 

UO 

V 

wc 

WGP 

WNP 


Definition 

Any  6  digit  number  used  in  the  random  number 
generator. 

Quantity  in  pipeline  to  be  delivered  in  K 
week,  product  2. 

Quantity  in  pipeline  to  be  delivered  in  K 
week,  product  1. 

Price  markdown. 

Uniformly  distributed  random  number  between 
0  and  1 . 

Salaries  for  personnel. 

Standard  deviation  of  demand. 

Stock  received. 

Sales  stimulation  foctors. 

Change  in  number  of  personnel  in  the  current 
week. 

Total  accounts  receivable. 

The  mean  weekly  demand  for  units. 

Total  units  sold  during  the  week. 

Units  of  lost  sales. 

Quantity  of  units  ordered. 

See  DEVIOT. 

Wholesole  cost. 

Weekly  gross  profit. 

Weekly  net  profit. 
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Variable 

WNPT 

X 

XM 

Y 

Z 

ZT 


Definition 

Weekly  net  profit  total  of  product  i  and  2. 

See  RND2. 

See  M. 

Number  of  salesmen  for  coming  week. 

Number  of  salesmen  for  current  week  product  1 . 

Total  number  of  salesmen  for  current  week 
product  1  and  2, 


DON  AFT  Subroutine 

The  DONAFT  subroutine  and  its  internal  subroutines  are  shown  in  the 
flow  diagram  in  Figure  4.3.  The  internal  subroutines  to  the  DONAFT  sub¬ 
routine  are: 

MULTR 

MEANS 

SUMSQ 

INV3  and  ADJUST 
MXMLT 

The  variable  names  used  in  these  subroutines  are  defined  as  follows: 


Varioble  Definition 

AAA  A  working  matrix  for  the  matrix  transformation 

into  intercorrelations. 


AUK 


An  array  that  contains  the  sum  of  squares 
of  the  deviations. 


B 


The  matrix  of  weights. 
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Variable 

BO 

C 

CUB 

DELPSl 

DELX 

DELY 

DET 

DTAPSI 

MR 

RSQ 

S 

SB 

SBO 

SEE 

T 

UB 
UK 
X6AR 
XX  (I,J) 


Definition 

Weight  of  the  Intercept. 

The  upper  left  quarter  of  the  CUB  matrix. 

The  inverse  sum  of  squares  matrix. 

Backward  difference  of  the  SlMCO  state 
variables. 

Backward  difference  of  the  SlMCO  input 
variables. 

Backward  difference  of  the  SlMCO  output 
variables. 

Determinant  of  intercorrelation  matrix. 

Forward  difference  of  the  SlMCO  state 
variables. 

The  squared  coefficient  of  multiple  correlation. 
Coefficient  of  multiple  correlation. 

"S"  matrix  of  the  DONA  model. 

The  weight  standard  deviation  matrix. 

Standard  deviation  of  the  intercept. 

Standard  error  of  the  estimate  orroy 
Student  t  statistic. 

The  upper  right  quarter  of  the  CUB  matrix. 

The  upper  right  quarter  of  the  AUK  motrix. 

The  meons  of  X  dato  motrix. 

An  array  of  the  forword  and  backward  difference 
of  SlMCO  variobles,  olso  X  data  matrix  used  as 
an  input  to  MULTR. 


DONA  Subroutine 


The  DONA  <iubrouHne  is  shown  in  the  flow  diagram  in  Figure  4.4. 
variable  names  used  in  this  subroutine  are  defined  as  follows: 

Voriable  Definition 

R  (1,1)  Vector  of  independent  variables  in  the  DONA 

model.  The  "1"  index  denotes  the  particular 
variable. 

RR  Same  as  R  (1,1)  except  the  "J”  index  is  the 

time  variable. 

S  "S"  matrix  of  the  DONA  model . 

T  (I  1)  Vector  of  dependent  voriables  in  the  DONA 

model.  The  "I"  index  denotes  the  particular 
variable. 

Tf  (I,1,J)  Same  as  T  (1,1)  except  the  "J"  index  is  the 
time  variable. 


*  FOh 

k.  pmcokam  FOH  ‘SImCo  *nu  OOua  cabohaTORY  System 

CALL  'jIMCO 
LALL  uONAFT 
Call  l»ona 
STUH 
ICiU 
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l  UAHS  I  l).tYll.4)).|FCSU).YYll.5)).  'SALS  I  I  I  )  *  YY  11.6)  )  >  IVCSl  11  )  >Y  f 
2U.7li.(CSSlU).YYU.B)).ltAMSlll).YYIl.9>).IFCSlll).YYIl,10)) 
EGolVALtNLl  iAr;0«SlU.Xll.l)).<Al>UiY5111).Xll>2l  ) 

{  oMk  Ol./HANO/lMD.  'NSTART  .HN01»  iPN.NOPS 
nST  AM  t : 1 V34io 

call  i’.n 

0  FOHMal  (///k>.4'i«t.rA.  4HAUUH » 4X  #  SHAOO*  l.bX.3MSAL.bX.2nwC.AX.ui'H;ii, 

•  o«  .  J"LO'..  S»  '  l"I'IS.6X  .  JhGIS.  7X.  iHVOP.  7X  .  5H»NP.  4*  ,  jM/jLS.  S«  .  2hFu. 
*.<m2  t  •  it  t  1.1.  .0X.2hmt.4X.Jm  S"/) 

L  l•MINI  COLUMN  MlAUINOS 
•  mUl  lb.  M 
L  (i.nlALUt  VALL'(S 

..  I'-IIIAlUl  UAIA  llCIIMulAT.ON  vALuCS  A1  2ER0 
..  ll.niXLUL  PlltLint  fOH  '.ME  *.txY  B  RtEKS 
DO  S  Mil.b 
S  PU'LlNlMlIiU.O 
Ul  iU  m: 1  .  b 
so  PlMLlLlM)=S(J.l.' 

P  IPL  l  .l  7  j  lO  . 

P1PL1L17I=0.u 

PlPLlNia)=o.u 


T».u  :  bO.O 

SuS^U.O 
buUiu.  0 
iJC.sbO.O 
UOl:bJ*n 
ULj  s  u.u 
ULbJ  •  J  «  (I 
AlL»-ibU.U 
Aitul^lbu.U 

Clw  "  t,u 

Ch\«l  su  •  0 
Lit-  r  ^  :  1  •  i  u 
o  ui-i'  <  J  i  >bu .  u 
J(  bU  1  • 1 U 
ri(J  I  J  )  •bo  •  u 

f  I.  -*>0 , 0 
►  (  )S‘JU«0 
u  -  0  «  U 

»*»  oii^y.o 
'  1*1^  •  0  «  <' 

'»‘i*  i  «  o  *  f' 

I  ACCiJuNtS  HtOllVAMut  THt  NLXt  tt  JfctE^S 

iji;  7  ^  *  I  I  M 
/  A»/(»ia<;voo*«j 
IJo  7u  K  S  1  .  u 

■  ■I  Ah  I  in  t  Z^'jhO  ,iJ 
TAk:i,jU0V4O 
r  AHUiOOuli.v; 

COH5'*jOO.U 

Cv.ni:sUU(>«0 

rilUui^L  «ri  A;.T|ONS 

bV'"**'  z  o«o 

Y ( l»sj,u 
Yi  (n:i.u 
ul:)ASo  *  u 
^LbAl  -*.0  *  u 
A!-  U  »  50  .  W 
A*^l  (  iiZO.U 
.fH,  fCH  Ivll  •tk.f.b 
IhSluo-l 

UC  eO  •  in 

Ca^L  hANLi^*  ( *' *  u  ) 

OALL  bO*KLY  70  utlfUMlr.e  LtAO  T I  rt£  FOH  UNIT^  OHUE»tr 
CAUL  iUHHLY(HircliV*bH«UO«H» 
i»CSSK*bO« 

MLL  CU^foMtH  llMOKOWUC^^b  If  ant 

i  Alim  « 

If  Ku\^  •  Aioml 
u  hb  5  ouu 

AAlli»  S  A]h»i  -  CUO 
Cbg  5  u.u 

~  li.U 

OC'  7v  iO 

'y  Hi,  Z  A  1 H  W 

AAih*  s  o.n 

Ut  J  -  Kliu  '  AlH«)  •  *20 
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CUW  :  (CbO  •  *1U«)  •  .ao 

10  N  :  I.sa 
C(<J  5  N 

n  s  i.lS 

UL!»  S  M 

^  l*lL  DkHAUO  TO  otTCRHlNC  OCMaNO  t'OR  THIS  HCCK 
CALI.  UtMANOIAlJQH>PMO>TMU,SO«SSF>*) 

V.  oLTCMHlNt  total  SALCS  OUHlKO  TMC  hCCK 

If  lAuOM  •  AAIRW) 

11  tOS  I  BS  •  AUUK 

Ijo  TO  M 

li  TOS  i  Ob  ♦  AAIH* 

V.  ukTCHktNL  BACAONOCR  OUCUE  LCAtSTH 
Ct<u  :  CHO  •  lADO*  -  AAIe«lt 

..LTCHi^lNt  OHObS  U.COMtiOBOSS  PHOUT.  AND  ENOIN6  INVENTORT 

IJ  Olb  I  100.  •  tl.o  -  PMO)  •  Tos 

Lb:u.b*0IS 

Af.O  i  ulb  -  bu.O  «  TuS 
Alt»  :  AA1B«  •  AQUA 
If  (Mltml  Iviliilj 
1<4  A]  t  A  ^  0.0 
13  CO'Ull.UE 

XIw:AIEa*SO. 
lo  EC:7T7.Q 

L  uETEKl'llit  CAbM  ul.  HAND 

17  Cc.cLuM*  AH  (  1 1  »  .5«GIS-SH»?>0.-1  .2/32.  »  AbO.  aAIBN-I  .Ofl/52.  )  aTAR 
t  lOoAtt  AllOuNTS  receivable  and  DtTtftHlNE  TOTAL  A.«. 

Ov  10  0  f  1 ' 0 
1(1  Al.(j)  i  Am(j‘1I 
Ah (L I  J  .b  A  olS 

TAH  :  ANIll  A  AR<2t  A  ARIj)  A  AR<4I 
i.  uLTEHvlNl  aEtALt  i.E  T  PHOE 1 T 

anh  :  aoP  -  (.2/32.>AbO.  a  aiHa  -  (.0ft/b2.)  •  TAh 
11  tA..OA  -.To  A  FO)  IP. 21. 21 

H  Hv...tCAji  uL“Ai(0  fOH  fHt  CCiHll.G  »l£K  ANO  UETEmhInE  IE  A  PRICE  a-A'ikOOaia/ 
H  lb  NtLtbbAHT 

10  CAiL  A0HCAb(OuP<EU»A00»» 

if  iaiEa  a  HiPLiNllI  -  1.5  A  EO)  22a22>20 

VO  Pvtjru.L' 

bbEso.U 
oO  fJ  2} 

.1  call  aOHCAS  IoOP.EuaAOOB) 

v2  1'»W  A  O.U 
bSA  A  0.0 
CO'1I1.(OL 

L  call  OkL/Lh  to  CeIEhmInE  ThE  NunUtR  OF  UNITS  TO  HE  OhoEHEO 
C*Lt  OHOLH  IHlPLl...oO.TuS.AltiA»EOAASB>OI.DSB) 

CALL  hAnUMIA.bI 

CALL  boPHLTlPIPLltASHl.OOl.a* 

■ClIbHlAbo. 

AlBaUAltalASHl 
|A iCuwl-AlUal I  SO. BO. 90 
.0  UbiSLOOl 

AAIbalAAlBkl-CROl 
CLUIAO.U 
ULbl :o . 0 
oc  Tj  lOo 
90  HbljAlBal 
AAlBalSu.U 


! 
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UL!>|:tCbOl-AtU«ll«.iO 

CbiJi:lCbei**lBkll*.80 

too 

CB<ai:(< 

N=WLbt 

OLIilSN 

CALL  klCHAND  tAOD«ll>PMOl>TMUliS0l>SSPl>A) 

If  (AUO«l*AAtUM>  )  llOillOil^ll 
IIU  ToSi:uSl«AOU«t 
00  to  liO 

liU  TuLtsb&t«AAlb«t 

CUol:CU01*(ALiUHl-AAlB«tl  I 
UO  OlSlsl00.*ll>0-PMOll*TUSl 
C&l:U.B*OlSl 
liGPi:oItil'5U>U*TuSl 
Att«l:AA|HMl-AOD<«t 
lA  lAtCkll  14U>l50ilS0 
140  AUklsO.U 
tbO  LONtlNue 

(1V1:aIE*1*^o, 

IbO  FCti?/7.U 

170  CO«i:oOHi»AKlU  )•  ,!)»0ISl-SRl*50.-(.a/5*.  )»iO.*AIb»l-(  .OA/'ii.  )  •TAMl 
UO  too  j:liA 
loo  AHl I JISAHl tO«l I 
AM(4IS.b*Glt>l 

TAMISAHI < 1 I*AK1 (2) «AR1 ( si tAHi 14) 

«^H|:•OPl-<  1*50,  (A ,08/52.  )*TAR1 

lA  IAuUkl-.70<A01l  190,210>210 
too  Call  aurcas(uupi>aui>add)i1) 

lA  <Al£«l«PlPLICIl)-l.b*ADl)  i20>220t2OO 
2(10  PPUISO.O 
b&AlSo.O 
00  TO  230 

210  CALL  AORLAS(uUPl,Fol*ADOtil) 

220  PPOlAO.O 
bSFlAO.O 
230  CONTINUE 

call  PEN&NL  (FO>A01<Y>T1«MT>2>21«2T>SAL,SAL1,SUMMiSUNM1 • JwkS.m.mI I 
uf'UATE  COm,ULS>«nM  WITH  PEK&ONNCi.  CONUISEPAT  IONS 
C(jMsCOM-A  C-bAL 
COMlSLOHl-FLI-SALl 
wNPTSI*NM*wnP1 )>AC>AC1-SAL>SAL1 

Call  ONULN(PtPLIE,U01«TUSl>AlCRl,A01,ASB1.0Il>nSBtl 

JWKS2:J«K<>«1 
t  PmINT  HEbULTS 

•PITL  (b>24l  JwKS2,AOOn1<«C,lO(4.(i1S.«NPT.AO 
24  AOHMAI</|T>eAiFa,2«9X>F<9,2,llA,F12.2,6A,A10.2. 1  OX  .F  10 . 2  •62  •  Ffc .  o  I 
WHITE  (b,3M  AODw,SAL.TUSl,ToS>WOP.ULS>2Ti2>MT>SH 
3b  FONMAI  (TX.A6.0>ax>F9.2»9X>A11.2>12x>F6.0> IOX>F10.2,10X.F6. l,bA,F 
>b.2,2A>F3.1,2X,I3,2X.FS.2) 
wNPTbt JwKSlrwNPT 
LIObl(i;:bO.O 
UOb( 11:50.0 
SALS(l):bOO.U 
wCb(l):3000.0 
TosSdCbO.O 
AOUwbd  )  =  50.0 
AOUWbllDSLU.O 
ULbSd):u.o 


UL^iUll:0.0 
CbSSl 1 ) £U.O 

CbU^Kt  i:0.u 

FbbI 1 I3b0,0 
FO!>1(1I£!i0.u 
rAHbui^tooou.o 
tah&i 11 )£toauo.o 

COh>|I  1 1  :«aou>u 

COMSl ( 1 |:40uu.0 
<&lll:0.0 
2&i(ii=o.a 
MTidlso.O 
bbiu  iso.o 
bb^l 1 1 1 :u.o 
AILi»bl  1  i:OfO 
AItlir<il  (  1  ):0.u 
Ubbll i:0.0 
MSblUI^U.O 

(11:0.0 
Aul(l>S(  1  ICO.U 
lUbSl ( 1 ) :0.u 
(>t!>&(ll=0.0 
OtbSKl  l:o.u 
■  1 1  :0.0 

•  OHii  I  I  )JC.., 

foiiusu.u 
FC  bl 1 1 1 :u . 0 

•  '.HS(  1 )  su.O 

•  f-**!!!  ( 1  I  lO.u 
ASoS( I ) 10. 0 
AbUbl ( 1 ) sd.o 
UIS( llsu.d 
(.<1^1(11:0.0 
USbSdUo.O 
OSoSl ( 1 1 :4.o 
Cbbdco.O 
CSm  ( 1 1  :u  ,  0 
kC  b I ( 1 1 :o . 0 
Salsi 1 1  I :o.u 
»  I  VS( 1 1 lO.O 

« 1  vsi 1 1  I :o .  u 
/tb(  1  CO.  0 
«;>(  ICO.O 
“SI  (  IcO.o 
OOS ( u«KS* 1 1 :ou 
UQSl  ( JaKS«l  CuOl 
utSSI  j»aS«1  Cues 
ULSSl  ( J»AS«1  CULSl 
AltOS(  OURS*  1  CA  le» 

AIti«Sl(J«ASdl:AIC»l 
Cfcosi-OAS*!  ccao 
CbuSl (0»AS» 1  I :C8o 
F I  S  ( w«AS  »  1  (  (j 
FCISI  io»AS«  1  cf  01 
lAMSl  J»AS*  1  CIAH 
lAHSl  i./»AS«l  I  :TAH1 
Ct.'HS(  J»KS»  1  I  :tOH 
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COMSl  ( JdnS^KsCOMl 
2S( )=^ 

2&1  <  jiiiK&«i  1:^1 

MTSI j«KS*l ISMT 

St<b(  Hit) 

SHSJ  <  )  sum 

JaiKS«i  IsAIBm 
AlbkSl lu«KS*l 

USSlJiiKS<i 
U&^l (WHKbtl ISbSl 

AUJB«ilvl«KS»l )  SAAIBHI 
AOOmS*  JNKS«U:«00* 

AUUh'jl  Ik^kKSti  ISAOOkl 
ToSS»jl«K$*l)iTUS 
IIJiSJ(JwKS*l)5TUSl 
6I^Si JkKb*l I=0IS 
CIS^l  ( JaKS*U:OtSl 

■  J»Kb«l  IStfBP 

•  Gf’Sl  t  .IkSPl 

FCblu«Kb«l ISfL 
PtSi  t  Cl 

■  SPi(  Juki*  I )  2|iNP 
kNPSl  I  JkAtitiiSkMPi 

Abd&(JiiKB*LI:ASA 
AbbSl • J»kb«i IZXSUI 
Olblj.kSM  (Sul 
01 bl IjKKb* 1  I sg  1 1 
Ubd&l  J«Kb<l  CUSH 

USdSl ( JkKS*! ( :OSbl 

CSbtj«KS»lC(.s 

CSbUuddS*l)sCSl 

■  CSIu<iKS*l  ISkt 

•  Cbl  ( Ji««S4i  c,ci 

SALStukKb*!  Cb*L 

SAtSUUkAi*!  CSAlI 
AlVSlJ«Kb«lc*IV 
AlVSi(JdAS«lC*tVl 
ilbl  J«AS*1  CJT 
“StukAS*! I sp 
Pbl  ( JkKSft  cm; 
db  COWtiKUe 
ML  ruHI. 

SobHOuT 1 Nt  utMANu I AOPk . Pwo • ThU • SO . SSF . V ) 
0  s  IHU  •  bu  •  V 
ir  «ul  1>2.^ 

1  0  s  U,0 

i  If  (P«0  -  o.lbl  3.H.4 
I  SbF  s  0.0 
oO  To  b 
■*  bSf  s  O.ub 
b  AUUk  5  D  •  U.O  ♦  bSFI 
N  =  AuOH 
A'JUk  :  N 
HCTuMiy 

SuUBOuTInE  b'JPPLT  (PlPLlN.SH.UOct 
OI**t'4blON  P1PLIN(9) 

UtTlPPlui.  SMlPPtldS  HtCtIVCD 
Sp  I  PlPLlNllI 
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t  uPUATt  buPPLY 
OCI  2  ><.  =  1 1  7 
/  PIPL1..IK)  : 

P  I  PH  111  a  1  :  u.O 
t  UtTtPMJNt  Lt*C  TIME 


MS 

XtiOO. 

If 

i 

LI 

S  4 

GO 

TO  iJ 

if 

S.b»6 

b 

LI 

=  s 

GO 

To  1  1 

b 

U 

(M-74)  7. 7. a 

7 

Lf 

=  o 

GO 

To  11 

a 

U 

IM-ttG)  9»V»10 

s 

LT 

=  7 

G(' 

TO  1  1 

0 

C  T 

«  0 

■  ‘tExS  Supply  OhuChED  BY  uMItS  OODEHEO 

U  PiPLl.ii-t  1  ;  PlPLimti)  «  uo 

W  t  f  U*’  • 

^\>«wooT  ir^L  ►uHCASIUOP«Fo»AOO*) 

uMt.iilON  LuP<  X  I  I 

E.U*^LATt  ►  00  Tut  LAST  T£i,  ft£iKS 

OL  1  lU-1*^ 

L  S  I i  -  Ih 

1  uc/P  ( u  *  ’  au*’- if  j » 

•yLt'i  l  )  _  A{»u*i 

V  LALCt'LAlt  ♦uMt'i.ASTCO  UtMANO 

►  wi  i  I  ♦uuPi2)  ♦i>OP<  5)  )/3,u 

HL  1  UH:« 

Si.'Ok0i.T  I  Nt  LHU£R  IP  IPL  1  Nn'O*  ToS*  AJtW/K' »  ASM  *  03  »  DSh  ) 

I.  *.AlCllAIl  OLsl^tC  iNvtUTOMY 
Cl  :  ^,0  *  ^  u 

V  C*iLCULATt  SuPPL  »  UACKQHOtM 

ASn  i  0.0 
lk.  \  .  r  I  .0 

J  ASo  S  ASM  ♦  ►'iPol.MN) 
w  '.•itCoLATt  ubSiHio  Supply  MAC*<o«Ot.P 
LSo  -  o.u  •  ►  ui 

Uo  =  /us  <  .b  ♦  iUl  "  AIt»  ♦  OSB  •  ASH) 

ujoo 

UU«M 

if  <  vv  »  f-  o  •  J 
i  UO  i  u.O 
3  Co.4l  i.>uf 
wt  fuM^ 

SudwuoU«-(L  «A/4C«IULVlOTfWND^» 

Co»*MOl,/MAf4C/  IHO»NSTART.Rnoi  ,  iMiM,^0PS 

I  WO;  HHN-l  1"N/I0  00000000l*l0unuu0000l  /luoo 
I.  1  .Lv«.  0  I  iMOr  1 

'.OP‘j:uOP',i  I 

NOOi^f  LOA  I  Unut  •Q.iJOOOOl 

Ot  Y1  W  I  ;bwKl  ,  .  •AtOG(R»iOl  ,  lAtOb  lo.Jfiilebn.RNOi) 

wNo  I  ;"ii.o<’ 

Mt  /  UiM, 


THIS  hElK 
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il>OHULT  INC  iMT 

CgMMCi,/HANO/IMO*NSTAf)T«RNDl»  iRN«NOP$ 

ImuSN'sTAHT 

HNv^l=*-LOAT  (  NbT  AW  T  1*0. 0000  0  J 

f*OPS  =  47'*Jti 
wt  rijHN 

SUbWOuTlNE  fr-LHSNL  «F0tF’0ltt*Tl«MT.7.^1»2T#SAL.SALl»SuMW,SOMMi 
•  •  H  «  M  1  ) 

uimln'jIon  \ '  toll  •  n  1 101 )  «iiH< loui  <xMi  ( 100) t kmt t  looi 

rijHK^«i  I  zfij/iO, 

T1  IJl«K.S*l  lst-Ul/20, 
u;=r  (wiiKS«l  I  *0.9094 
Li:YljiiK!il*0.9949 

uUzrKJaiKS*!  1*0.9994 
Lit  =  ri ( jkK^i*Q,4qq9 
Ml=H«*Lll 
«*'  I  JUKS*  1 1  iw 

»*.t  (  J«Ki*  )  I  :mi 
SOM>'-bU“*'*X|..  I  JKXi  I 
■jl'M**!  =bi)HMl  •  »wl  IJ«KS  I 

/:riuMM»r  1  1  I 

<;i-t.uM*ii*Tm) 

HTSMl *K 

.•b  «*  I  I  J«*.b)  IMI 

It  (  »..( j»*bi  I  l.J.i 
1  bui.-;*:  •  19«! . -»**ij«Kb  I  •* .  *19.; . 

00  To  0 

i  b«L:<l*l0i. 

«  It  I  XMl  I  J«Kb  I  I  b.0.6 
5  bAi.lIil*l92,-«m(J»(HSI«?.«19x. 

00  To  7 
o  iAn:tl*lV2. 

7  lit  I  UH. 
ti  0 

.  (  OH  UONAt-  I 

bUOWuoTINL  OUHAFT 

.OAMAMt  (tt<  I  »»«:2  .  I  ,  Ifsi  iJo.  10H:ft2.  IfiOtP;  7  .  ya«t  lo 

C0itA'oi./8L«  1/(.UH(MA«,MAU 

Co.tnO./tlUK^/PSl  (  10e>  IPSI  I  ,  »  I  tOdr  1*X«  I  »  YV  1  lOui  I  TY»  i 
CO«MOU/Hl.Ai/il  KtOtP.  INCEP) 

OIMCKiilON  OttPSI  I  lOB*  IPSn  .01  APSlUOe*IPSJ  l  .OElXIIOR.  l  YTYI  .OEL  r  I  lo 

lo.lTtri 

UlMENblOO  »t<MAK.Iue' 


I.  r.'L  PULLO.ltjO  calculates  del  pSI. delta  Pbl.UtL  A.DEL  » 
IbS|Oo-l 

00  lUuO  1-1 . IPSI 

00  |OiO  lU 

OtLt'Si  (  J.  I  l=PSl  (  J.  1  I -PS  I  1  J-1  ,  I  I 
OTAPbi ( J. I  lIPbl (0*1.1  l-PSl ;o. I  I 
lUlO  CCNTlliUE 
lOUO  CONTI. lUt 

OO  1UU2  l=l.iAXX 

00  1012  j:2. ID 

OILX (o. 1 1 :x ( J. I !-i (0-1 . 1 ) 

10)2  CCNTlNOE 
1U02  CUNTli.Ut 


99 


UO  lUOj  UlillTY 
DO  1013  j:2>IU 
OtLXji  1  )=»MJ.ll-YY(J-l.l) 
lulj  CONYlNuE 
iuU3  continue 

■HITE  lb>20JU) 

20CU  FOHHtT  I  IHl  .tiOY  f  '  T  I  ’  /  13X  i  IH?  >  61  ■  IH  j  i  eX  •  IHH  >  BX  .  1H$ .  9*  •  1 H6  •  SY  i  In  7  , 
laxi  Inti,  ax  .  Ia9<  7X>2Hl0iTX,2Hll.TXt2H12>7Xt2H13i//.3SX>  OEL  E-SI 

2.  aALMiAHC  HFfCSENCE  OF  SlAIE  VARIABLES  — — •) 
mRITe  <t>i20ul  I  (OELPSl  ( J>19)  •u:2<l3) 

^uul  FOMMXI  (•  TuS* . 3ltl2F9.2» 

■  RITE  I  6  >201.2  I  (DElPSI  (  J>9I  .UAa>  13) 
iOu2  Fuhmat  <•  fl')9X.12F9.2I 

•RITE  (b>2003l 

<-uul  FoHxAi  (///41A.’ - DEL  X.  backward  DIFFERENCE  OF  INPUTS - ') 

•  RlTt  (c>>20uK>  I  <U£lx(J>k  I  ij:£>13I>k=i.2i 

2UV>>  FoMMAI  (•  AluwS*  t  U>  12F9.2./4  •  AOUWSl  •  1 12F9.2) 

•RITE  (6>20uSI 

FuHwAi  (///jax.* — —  delta  RSl>  Forward  niFFERENCE  of  state  vamia 

IMLES - ’  I 

•Rlrt  lo>20olJ I0;aPSI(U>19I>u-2>13i 
•RlTt  I6>2Uo2I<OTAPS1<U>nI >U^2> 13) 

■HITE  (b>20uDl 

.,.1.  fuRRAi  i//,i,iA.> - OEl  t.  backward  difference  of  Outputs - - 

1 1 

■  HltL  lo>20i,/)(IUFLTlJ.KI..iii.l3l>Kil,2) 

,  1.  7  FuHMAl  (•  sal  ■  >  3«  .  12F9.2./ .  •  •£>#12F10.2l 

L  I'l  FlllLaI'.U  FlKMS  r>.t  X  data  wkTkIx  Fhuh  Sl^'CO  DELS  AND  PfllAS. 
uf.i  1  uu-2 
UL  lJud  UTl.M. 

XM  (  1  I v)  AuLLHSl  ( J»l .  1  I 
AX  (2.vliULLi-Sl  <  j‘l.2l 
XX  (  3io)  lULLI-Sl  (  J*l  >SI 
I  1  I L . J I :u(  LHSl  ( J* 1  .H I 
AX  (  '  .  j|  Suf.LHSI  (u*l .  10  ) 

X  < I e • U I :UELA ( J« 1 • 1 1 
AX(  /.j):utLA(j*l.2) 

AX (  o> Jl :urAHiI ( J»i . I ) 

Ill  -I .  u  I  vuTAi-S  I  (  U*  1 .2  • 

X  X  I  1  I.  .  J  I  Xv,  I«HS  I  I  ^*1  <S> 

XX  (  u  >0)  -•..TArSi  <U«1  >9) 

AX  (  Ik  .  u  I  LUI  At  ;  I  (  J«1  .  10  7 
AX  U3>  Jl  Xl.t^T  I  J>  I  .  1  I 
A  X  I  I  L  .  o  I  XUl.L  I  lu«  I  .<■  I 
.  C  0'>  T  I :  .UE 

*H  I  r  t  <h . ROi i 

oul  FLHHAt  !i)il.SSXf>  I  HATPIA'I 
1d-“Ai 
JU=N''' 

Ku:  (  1  iu-  1  I  /i‘ .  •  i 

■  HiTL  (b»  1Cj)<  (X*  J  l»J>  t  .J3i  t  J\t] 

i  :  1  f  OKfc'A  i  i  i  «  » fly  lb  ,  /  J 

•  i  Tl  U)  •  1  r  I  ) 

I  t  i  ►  CMVA  !  <  1  FT  I  I 

Call  '•ul  t  t  a  a  #  i  Ni;t  x  1  nOL  a  i 
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(.  Tut  fOLLutilNO  T*HES  THt  RIGHT  MALE  OF  TME  "CUB"  MATRIX  AND  TRANPObEb 
L  III  OlviiiG  tiiE  "S"  matrix.  TtiE  "S"  MATRIX  IS  SIZE  INOEP  X  INUtP 
UO  UIJ  UI.lNOCP 
ir.Tzli.OEH*! 

UU  j:lST.MAX 

iHi  SIJ-l>.OEH.Ii:CunU.J) 

■HITL  16.44441 

44U4  F0HM41  (IHl.bSX.'  S  MATRIX" 

iiHITe.  I6.444S)  (  (S(  I  ..j)  >u=l  •  INUEPI  •  1:1.  INOEP) 

HsuS  RC'HmAI  Um  .7G)3.SI 
ML  I  UHi. 

LNu 

.FOR  mijLTM 

SUbROuTINL  POLTM  ( « .KK4 .M.N.MAX ) 

F'AMARLTtK  KoLS  =  aU.  IRORS:  14  .  l  =  T 
HLAc  *4 

I  UPU»i)  MUST  i.OT  b£  greater  TmAn  (ICOLS)  FMuiO 

OIMInsION  C iLiL ) >UK (L.L I .uH(L iL) 

ClMmu,,  a  (l.L  I  .  AR  I  IB0«S.  IRORS)  »  AHSO)  IROkS.  IR0«5)  .XBAH  (  lR(j» 

IS) • auk ( 1N0«S> IR  TaS) .see ( IROmS ) .SUO ( iRCaS) >  bui I 

XRU«SI  •  Tu  I  1RU4S  )  .  MR  ITROViS) 

UlMft. ilON  Sbi  1H0«S>  IROwSi  >  T  ( IRUaS.  IROwS) .  AAAI  IROaS.  IROwS) 

01  mens  I  ON  X  I  iROvS.  KOLS)  iRI  IbOkSi  IROaS) 
t0MMu../HLK  I/CUBI  IHOaS.  iPOaS) 
tool  valence  lAiB) .  1 Tfl.sei ,  1 AhSo.T ) 


^  IL)  :  '.0M)1fR  independent  vAHIAoLES.  (M)  :  NUMBER  OEPtNOEuT  k.kuIi, 

u  VAHIAuLES.  IN)  r  number  OOSEhVAT  ;.0NS  '.■MuTii 

If  il.lt.ni  oo  To  no  .•■rusu 

■  I'  J  1 1  I  6  •  IZ  0  I  1-  R  1  Uij 

tuu  EoHMA  r  I  TS  .  ’  tRROH  -*  NOT  ENOUuH  DATA"  ‘t-jlu 

S4T  stop  ‘'RsSE 

1  10  Coni  |..Ot  ''R.  ii, 

u  DATA  ITEMS  I'EAri  IN  ONE  OBSERVATION  OF  All  VAWIAHi.LS  AT  *  ’RiS. 

..  time.  <*)  IS  THE  UATA  MAIRI*.  ''laai. 

LAuL  'tAfiS  IX.  XbAH.  MAX,  N.lHoab) 

call  ,>umSo  I  X.  Aiik  >m/ x.n.  IRuaS )  '  RlaO 

IXbAm  contains  The  MEANS  OF  (MAX)  VARIAHlfS  (  AiiX  )  C  ON'T  A  I  Ns  "RJni. 

L  The  3uMS  of  SooARES  OF  THE  OlOIAHONS.  lx)  CONTAINS  InE  ‘■mIn. 

C  UEvl-MONS.  (SEE)  IS  IHI  STANDAhO  ERROR  OF  ESTImATI.  (CoH)  Is  i  •'!  ••RlSv 

L  INvtRSt  SN  matrix.  (B)  is  Tml  MAIrIX  of  aEIGHTS.  ISB  IS  T.<E  «£1»..T  'RiBO 

C  STU  OLV  matrix  'Rlv.' 


aRITl  Ib'VlSI 

SIS  FUHMAI  llHl.euX.*  AUX  HAIHlX" 

JOl IMil 

lUl^MAX 

KulSi  I  001-1  l/l))*l 
00  91J  KIl.KOl 
*H:a«N-7 
APirbxK 

aPlTL  ifa.lOvjiuAoml.Jl.j^KLi.KRn.lil.Iul) 
aP  i  T  E  1 6 .  1  C 1 1 
Vl3  ConTIi.oE 

DO  4l.S  Ul.RAX 
JO  41G  j:l .pAx 
4  1  i  AAA  I  1  .  Jl AAUn I  I  .  J) 

L  I.Kt  foLLO.lNO  REPLACES  ZEROES  allH  EPSILON  ON  The  MAIN  DIAGONAL  Of  ajk 
IF  I  AAA  I  1 . 1  I  Ibl3. 8|4 .R14 
blJ  IF  ( AAAI  1 , I  I .UT.-I ,E-6)  AAAIl.  I):-l,E-b 
00  To  ,130 
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"^3  *•*<!. I, =l.e.6 

■KHl  (b.<Jlo» 

■<lo  (hOXi>  ttA  MATRIXM 

CIO  Vi-,  xzi.nijl 
•'L2i()»R-  f 

llx  I  VUI 

^  «N,tRC0RHf  *noN  coerFiccNrs  o.., 

t*'!  I'ju 

^  "  ‘l'o  -l.£-6.5£T  IT  eou*L  TO  l.t-e. 

CO  iUu  l:i.>«A» 
j'  lOoBl  1 , 01  I  All, BIO, 

•>>  I  It  (h.hOl  I 

■H.I  K«m*i  <U.l,Ah,.MNttHCCF«luRTI0I,  «AIR1,., 

I  k,£M4  Jl 

*  VS  <  1  JU-  I  I  /<•  >  •  1 
w  V  1  u  *  K  :  !  ♦  K  Li 
7 

i  u  (£  .1.11) 

i  ''  I  ►  'r  H.S  «  1  (  J  ».l  > 

vV  M^O  1  •  1  «  L 
bcO  vi  I  .  L 

<'.  (•  *11.01  sCi-Ti  I  I  >  o  I 

o  All  1  (,v  j  I  ■' ' ^  . I  . .5(.  1 1 

•■'  I  '0  (b.l-.lTv  l.,f  r 

. .  ‘-V£R.t  MATFJ.  *  MSTI„CO,.aArn,  . . 

••'irc  io..-no...(.Cil.JI.Jit,i.,.u,,,  , 

o'o  fwn.Ri  (  i„  .ioi3,5i  ‘  *■  ’ 

w  dooM  4  ;  1  . L 

t)oo6 

,„VL«S.  HATPIX  *  1,0  .NTr.cappE0*,.0„s, 

*0110  lO'oSl-JI  (  (0  I  I.o)  ,0zj,7j  ,  Jjl,7i 
111,1  HHyAI  u„  ,7013.^, 

oc  jujiil  l:j,t 

00  lUyOl 

lO'JI  UM  I  ,  „  I  ,  , 

t.i.iv  FON..4I  ll..l,,!o*.,  R  IR*tH|,.) 
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xHtTE  («it6au^)  ((LIKII>J)iJ=liL)>lSliLt 
CALL  MXMlT  ICiUKiUe<l.>liLtL>L> 

MHITE  I6i6eibl 

oHib  AOHHAI  (IMltiiaXi'  UB  M*TRJX*» 

wMTL  (b>6auJ)  I  (Lb  1 1  ■'^1  *  1*1 'D 

HlCONbTSUCT  THE  Cub  matwix 
Ol  IAIil 

UO  j:l*L 

o<:<i  CuU(tiJi:C(l>J) 

00  eio  iii*L 
CC  Bio  JoliL 
Ulo  Cobl 1 • J*L I=Od ( I • J) 
tel  NS»:C) 


IE  ((,.t(3.LM  i  «'aj»o 

0  CALCOlATINO  bTAhOAKO  C'EVIATluNS  AND  INTEmCEPTS  ‘WJoO 

IMT=l*1  ►RUBS 

UO  iiij  U=lHniMAI  kRJ70 

K:.i-L  wRJ7b 

UO(A)  I  XUAt<<Jl  ••‘Rj7fl 

SLE(Miaua  (jiUI  vhjbu 

00  ilu  lAl.L  MHJOO 

UU(K):tJO<E  l-LuU<I>JI  •1BAH(  I  )  E'RlRb 

<■10  bCt  l»  l=SE£  (A  l-Cubt  1 .  Jt»AOA(  I  .U)  •'M-iOu 

HR (X isl .-(SEE (K 1 •(M-l 1 1/(AU« (U»J)*(N-U-1 1 >  ►Bo  Ob 

bit  (M=buM7  (AhS(SEE(x  l/(»«-L-l»  1 1  iHlIO 

blo(M:btt(xl/SOHT4(4l  **BLlb 

T(j  (A  I  surj  m  I /aaO  ( A  i  mRoIE. 

oO  tAu  |:l.L 

B ( ( .A  I :cub ( 1 • J 1  ' Po^b 

ill  l  I  !►  I  :bt£  (A  I  AbUHT  (  Aab(Cutt(  i  •  1) »  » 

til]  T I  1  lA  1  SM  I  I  >A  )/Se(  I  .A  I  OfJA 

L  STAim  of  E'BInTOuT  «mtu 

tIO  00  3uu  UFl'A 

BhU.T  ilO  vBtbO 

BBiNI  Hu 
PBIM  HO 
Jlo  FoHPAIUKOI 

IE  (iJbA4tO,OI  60  TO  J20  iKobu 

PBIM  i3o>U.uo  I  Jl  •  I  I  «B(  I  (Ul » isl  *0  )  l■'^ll.U0 

100  FOHMA  (Tii'blNCE  THE  UIjMHER  uF  vARIAbOtb  JlibT  EoyAOS  TME  I'lJuntB  Of  f-RsOO 

1  l,ATA<  NO  VAHIANCES  can  BE  CUMPUlEO .  • // T  i  .  •  b  (  R’-Tli-LlL  fWsuO 

i.a/(  li.  *U1  •  .  li.  •  .  U.  •  f  S- .  I  li.tlB.b)  I  "RbOO 

oo  to  300  EftblU 

JtO  BbUAboHT  (  AbblMB  (U)  I  »  AMbU 

PBiM  3Ao*J*BN<ui*HSa.srt(jt#po(j».stiO(ji*To(j]  .ii.uii.oi.b 

lb(  I  »o»  <  T(  I  .w)  .liI»L»  AiBaiu 

jwy  EOHMAf  (T,i. 'OtPCNOENi  VAR  1  A0Lt  • «  T^Of 'MULT  I PLE  H-SOD  EbT  •  •  •  T40  •  •  Molt  f  Bblu 
iTlPLt  H  tST . • . T60. 'STO.  ERROm  OF  tST . • » /To . 12. TJO . 3 ( 614 . a . fcX I /// RHaiu 
l-'iHOlP.  VAKlABLEHTJO.'REIOMT'.rvO.'STO.  OEV.  OF  hT  .  •  >  TbO  •  *  blUOtN  mBjIu 
lT  T'i  rv.  •  lNTtBCCPT'.ViO.HOm.$.aX)/(TR.12.T20.S(Oli..a.6*  > )  I  MHb3 

300  CC'Nri'.Ot  fRboti 

Mt  TuN]. 

INO 

t  Fob  MiANb 

bUbBOoTK.e  BtANSlX.T.MAX.N.jBOhb)  Ml, bulb 

t  COMPorti  TMt  MEANS  OF  IMAXIVARIABlES  RAsiO  ON  IN)  OUSE  BOAT  lOfiS  «(■!  M-.Soi" 

0  Ett  ANO  rIEoOS  DEVIATIONS  FROM  THE  MEANS  MN'.jtl 

OIMENSlON  A  1  IH  3wb.NI  .T(  iBOAb) 

00  lUU  1  =  1  «  max  M.iSix,- 
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« I  1  I 

OC  1  lu  •!< 

IK  >  I  1  )::ri  1) 

Y(  1  t:i  (  n/N 
uO  1  ub  k:  1  > i< 
luu  X  1 1  .n  I  :i  I  I  iK  t'T  •.  U 
«t  T'.ihh 
tNU 

Y  bH  SUHSu 

^UHHuuT  |NL  iiUMSQ  (  >  i/tiMAX  <N«  IxOd&t 
FIHLi  THt  SuH  0^  ibU*«ES  Of  l*J»  OUSEftVAT  tOMS  OF 
tyilh  in  a  ^cpanate  mathix 

DIME  Nil  ION  X(  IHOa&iNl  •A<  INOWSiMAXI 
UO  lUu  KlfMAI 
UO  lUu  J=l I  1 
A( 1 • Ji:u. 

DC  tiu  k:1  iN 

ilo  Ai  i.jI:a<  1  I  jl<IU  <KI*XCJiaI 
loO  Aiu, 1 j:A( I Ibl 
Bt.  ruwN 
C'lU 

Y bo  MlMLt 

SilliBUo  t  /  nt  0-x«Lr(AiM«C»M»L»NiNXiLXJ 


MNSObO 

MNSUbO 

HNSuTU 

MNS060 

MNSU90 

HNSlOO 

MNSllO 


■jMbl/UlU 

IMAXIVABlABLEii  AND  i>M<>0u2u 
iMsou;.'_ 

ijMbOOAA 

sBNOObO 

ilMbOUbU 

bM&OUTO 

bMbOUtiu 

iMSOmo 

iiHiiOlOO 

MBsbO 


BtAc  xAIUl.  MUL r IPL IC at tON. 


UlMLN  jlUI<  A  (MX  iL  I  •  tt(LX  iM)  .1  (MXiN) 
UC  1  bSl.f. 

UO  1  Ul.« 

C  (  1  .uKU. 

00-  1  OII.L 

1  c  (  1  .  w)  :c  i  1  .  ul  •*  II  .  Jl 

ot 

LOU 

1  \  o  (  f .  V  3 

ioaUouIlNt  i N< K d • A iN.OE I  I 


oI'OhOE  TOWNSEND 
UEPT  I'To/NDO.  6H9!) 

S^dwUuTlNt  Yum  iNvtuTiTiO  SO  ,-«DE  MATRICES  wHICh  ABE  N  HT  ft  OR  LtOO 

I  to  i..¥EBT  l.AK(,tH  MATT- ICI  S«S«T  M  dT  Ml  DIMENSION  IP  I  irfOT  I  o  I  .  A  l  M  WM  ) 
If  ut  »  I  >  n  iCT  i  M  )  ANO  NtCOMPlLt) 

I.  Is  <Nt  SUL  0^  A  I*  IS  an  or  N  mathix) 

Til  SuULOurU.t  HLTUMNS  A  INVtPSE  IN  PLACE  OP  A  AMD  THE  OETE«olN»l, 

If.  uLIEMA, 

01  ME  NS  I  ON  IPI vOT ( 7) 1  INDEX ( T.^t .PIVOT (7) . A  I  I ) .yl 1 1 

II  Is  INt  SUuAkE  math  lx  TO  be  INVEHTEO 
A  IS  iHt  AnjoSTEU  INVEHSE  in  at  n  i 

LoulvAOENOt  llBOn.JHOwl.  I  ICwLLiM.JCOeuP)  .  IAmAx.T.SWAP) 

A  - 

UC  1  1  A  I  '  “ 

1  A  I  I  I AU I  1  I 

f'lOl  AbAl  .  ... 

Ov  <;  I  A  I,/, 
i  U|U  1  Av,  =  p  U  I  Au«PI  I  I 

L'l  r Emma  1  .  If 

to  ^0  0-1  ri. 


luVHOO  V 
li.VHOOlU 
li.VWOUSC 
4  : 1  V  o  0  0  fl  IJ 
li  ,  ¥  A  0  U  7l. 
1  I  .  VU'JUrt'i 

4  ,  VH  O.i  f' 

1  .VHOI  lo 
r4-,voui.'.. 
ii.vMOioif 


I..V00100 

il,VR0l7(i 

liiVoOiSii 


i  ,VIm'.!1iL 

4  ..VI.  1  >1 
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<0  IPlVOUJIsO 


DO  SbU  1=1. N 

lNVHU<2b 

*MAK:u.O 

INVR023U 

00  lUb  JSI.I, 

1NVN0240 

IP ( IPIVOTI JI-1 1  b0.103i60 

I'NVPOibO 

oO  00  lOU  KSl.N 

l..VP02bO 

IP  UP1V0T<K1-1|  a0,100.7K0 

1NVH0270 

bO  M  =  N*(r-1),j 

iNVHO^ao 

as  lEoliJS***''**''*”'*"'"’  •***00*»0I> 

INVR02<1U 

INVROJOO 

tCOLUM=K 

liiVPOjlO 

AMAH  s  A(M| 

1I.VKDJ20 

luu  CONTIi.ue 

ll.VHOJJi' 

lub  continue 

INVmOJUL 

IPIVOI  (  lCOHjP)  =  IPIvOT(  ICOLUM)*! 

Invmu  ball 

I»  (  IMu»-lCOtUPt  IbO.ibO.liiO 

liNVNOJbO 

IHO  OtTtahi.0£TiHi<l 

1 xVHO j  7u 

DO  2U0  L=1<N 

li.VXOjOL. 

M  =  N.it-ll 

ibVhOAN. 

“1  =  K*J{.OLuM 

1  » V  •.  D ..  j  ' 

i  7  V  j  W  i  C 

S»Ai.  :  A  (»'  1 

A  .1  V  »A  .J  •  i  . 

AIM)  :  AIMl) 

1  .V^OhJv 

2UU  AINl)  =  ^wA.^ 

..ou  IKUEAI 1 . 1 )  :mow 

iibVhOHD. 

lNL)t«ll.2l:K0LUM 

i  1  •  V  j4bli 

M  =  N*  (  UOLua-l  1  ♦  KotuM 

iiiVK0%7^ 

PI VOI 1 1 >  :  .IMI 

i '  4V*M  Oho  u 

L)t  TCnmsulTEkmaPIuotII  1 

AIM)  ;  1  ,u 

*  1 .  /  p-  '  j  0  0 

UO  Jbu  LSl.N 

I.aVmO^IIj 

M  S  N«(t-l  XICOLUM 

ImVp-  Or><o 

JbO  AIM)  ;  A(M)/^'^vOt  (  )  1 

l'.v»<Ob.iU 

DO  bbu  Ll  =  l.ri 

U.V^aO  jMli 

IP  ai'icuLUM)  Hnu.bso.Moo 

l-.VKOb‘jU 

H'VO  M  =  |.•(KOLoM.l)♦^.l 

i'bVnUSfii' 

r  S  AIM) 

4'»v-'r)‘j7u 

AIM)  :  0. 

T'lVfpOjoD 

00  Iib>j  l:1iN 

«  S  !<•  IL«1 1 

*»VhOoIjLi 

MI  2  -MtOLUM 

inVKUolU 

M  2  H.Ll 

li.VHOoib 

HbU  AIM)  3  AIM)-AlMl |«r 

li^VKOoi'J 

bbO  continue 

IhvhOd^ii 

UO  711/  ln.i. 

i .  iVK  Uobu 

1.2N41-I 

1  «V»AOot>«J 

in  iNutxiL.D-lNQtxlL.J))  6Ju.7l0.fc3o 

OJU  JMO*.=lfiot»  IL.I  I  *u»a30 

i  I'j  Ju7u 

Ai.VpiOuuo 

UCOLUMSlNOExlL.a) 

iNVMOo^iJ 

M  =  N*IJMO«-ll 

*  1 4  V  K  u  7  iJ  0 

Ml  2  N•(JCOUUM-^) 

li,vk071'j 

UC  70b  K21.M 

iNVKQ  /^u 

M  =  M«1 

II.VK0730 

“1  2  Mill 

i.bVKO  /!<(! 

bhAM  =  AIM) 

iUWKO  /'Ji. 

AIM)  :  AIMl ) 

1 » V  p>  U  /  Ml 

*l«l-'  :  b»AU 

ir.vwo  77i- 

7Jb  CONTINUE 

4  *VJxU  fn.. 

Til'  Continue 

i  1 .  V»i  u  7*^1 
*  N  VPb  0  boi' 
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InVMOulO 


li  vHOuio 
1:. VI'  in'll; 

1-  ,1  'li’l-: 


l>  £  TOiUStMD 

NOPTH  AMfnICtN  nCiCKWML 

I. 
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